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ABSTRACT 

THE HEAT AND SALT BALANCES OF THE UPPER OCEAN BENEATH A 
SPATIALLY VARIABLE MELTING SEA ICE COVER 


Daniel Reiner Hayes 

Chairperson of the Supervisory Co mm ittee: Affiliate Professor James H. Morison 

School of Oceanography 

The aim of this study is to paint a picture of the evolution of the 
horizontally variable ice-ocean boundary layer throughout summer. Observations 
were made during the drifting Surface HEat Balance of the Arctic (SHEBA) 
experiment in the summer of 1998. The ice-ocean boundary layer near leads is 
studied with an Autonomous Underwater Vehicle (AUV) and a novel technique to 
use vehicle motion data to calculate turbulent vertical water velocity along the 
vehicle path. Vertical fluxes are obtained and extend from the energy-containing 
wavenumber range and continue into the inertial subrange. This study is the first 
to measure horizontal profiles of turbulent fluxes in the ice-ocean boundary layer. 
AUV data are used in conjunction with: fixed-mast turbulent fluxes at discrete 
levels in the boundary layer, Conductivity-Temperature-Depth data (vertical casts 
and lead surveys), and a suite of measurements made by other investigators. 

The results indicate that scalars and their fluxes, as well as vertical 
stability, varied in the horizontal. AUV run-averaged turbulent stress in the 
boundary layer agrees well with the free-drift estimate. In early summer, fluxes 
were weak as ice velocity was low, and fresh meltwater was trapped at the upper 



ice surface. Also, surface melt was focused into leads rather than entering the 
ocean uniformly, resulting in a highly stable fresh layer. Near the end of July, a 
storm flushed leads, and the mixed layer freshened and deepened. The AUV 
observed strong fluxes under and downstream of rough, ridged ice. After the 
storm, heat and salt fluxes were strongest under leads. 

The results are simulated with 1-D and 2-D time-varying numerical 
models. The 1-D model produces a shallow, overly fresh mixed layer during the 
storm period. Simulations from the 2-D model suggest mechanical forcing from 
ice topography and a dynamic instability near downstream lead edges may 
enhance vertical mixing. AUV data agree well with the 2-D model after the storm 
and suggest mechanical forcing is important. The timing and strength of 
meltwater flux and the horizontal variability in interfacial fluxes have 
implications on large scales through the seasonal cycle of mixed layer depth and 
the surface heat budget. 
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1. INTRODUCTION 

This study investigates the redistribution of thermal energy and fresh 
meltwater entering the planetary boundary layer of the Arctic Ocean in summer. 
The fluxes of heat and salt at the ocean surface vary greatly over distances as 
small as a few meters. While openings in sea ice (leads) typically only cover a 
few percent of the surface, they are responsible for much of the energy exchange 
of the ocean and atmosphere during winter [Badgley, 1966; Maykut, 1978]. 
Because ice forms in the leads during winter, they are also sites of the strongest 
brine rejection. In summer, leads are the “windows to the ocean” since they 
absorb over 90% of the incident solar radiation, compared to transmittances of 
less than 5% for bare ice surfaces (even less for snow-covered ice, more for 
ponded ice; see Light et al. [2003]). The manner in which fresh water from snow 
and ice melt enters the ocean is dependent on the ice permeability, which varies 
considerably over the summer and within a single floe [Eicken et ah, 2002]. The 
latter study also finds that ice permeability also affects the coverage of melt 
ponds, which in turn contributes to the variability of transmittance. 

Numerical models of the planetary boundary layer generally assume 
horizontal homogeneity, a poor assumption for ice-covered regions in light of the 
variability due to leads and different ice types. In large-scale Arctic and global 
models, the fine scale structure of the ice cover cannot be resolved, so 
parameterization of the effects of spatial variability is required. Parameterization 
is impossible without measurements and understanding. For example, a naive 
parameterization of surface salt budget during winter freeze up might be to 
assume that the saline water formed in leads is uniformly distributed throughout 
the mixed layer. However, wintertime lead convection has been observed to be 
highly variable, depending on ice speed and surface buoyancy (salt) flux [Morison 
et al, 1992]. In typical Arctic conditions, one may observe either a “free 
convection” or a “forced convection” regime. In the latter case, the mixed layer 
becomes more saline and remains well-mixed, with a potential to deepen the 
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mixed layer under strong surface stress. The uniform mixing parameterization is 
reasonable. Under the “free convection” regime in which ice speeds are low, 
saline plumes do not blend into the mixed layer uniformly, but sink to the mixed 
layer base, perhaps decreasing the mixed layer depth. The uniform mixing 
parameterization qualitatively fails. This phenomenon has been simulated with a 
two-dimensional numerical model [Smith and Morison, 1998] and is consistent 
with observations [Morison and McPhee, 1998]. It is a direct result of spatial 
variability of surface buoyancy flux: if the total salt rejection were distributed 
uniformly over the surface instead of focused in leads, it would be much weaker, 
and the criterion for convective instability would not be met. That is, the lead 
number of Morison et al. [1992] or the Rayleigh number modified for turbulent 
diffusivities as in Smith and Morison [1998] would be too small. The effect of a 
strong flux in a relatively small region is not the same as a weaker flux over a 
proportionately broader region. Without measurements of the horizontal 
variability in the mixed layer, it would be difficult to assert that wintertime lead 
convection results in a shallower, more stratified mixed layer. 

The effect of horizontally variable fluxes on the planetary boundary layer 
has not been studied in detail for the stable case. Horizontal variability of ice- 
ocean heat flux over a portion of a single floe (less than 1 km ) has been discussed 
by Wettlaufer [1991]. He attributes the variability to coupling of topographic 
features with the structure of the boundary layer turbulence. It has long been 
known that the ice underside is not uniform but consists of relatively smooth areas 
punctuated by ridges and rubble fields. It is, of course, known that these features 
affect the momentum transfer from the ice to the ocean via form drag [Steele et 
al, 1989], and that variation in ice roughness leads to widely variable estimates of 
Reynolds stress at a given depth [Morison et al, 1987]. While the variation in ice 
roughness is important in the overall ice-ocean momentum transfer and ice-ocean 
heat flux, the variability in boundary conditions due to open water during su mm er 
is even larger. As will be discussed, the effect of summer leads on the planetary 
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boundary layer is also relevant in large-scale and climate studies of the Arctic, 
and the measurements and models presented here are relevant to those studies. 

It is well established that solar radiation through leads in summer warms 
the polar ocean mixed layer [Maykut and Perovich, 1987; Ebert and Curry, 1993; 
Maykut and McPhee, 1995]. However, the exact fate of this thermal energy is 
essentially unknown. It appears that much of the energy often goes into melting 
of sea ice [Maykut and Perovich, 1987; Maykut and McPhee, 1995; Holland et 
al, 1997; Nihashi and Ohshima, 2001], but the proportion between basal and 
lateral melting of floes is unclear, or more generally, the distribution of melting 
among the thickness classes is unclear. The quantity and residence time of 
thermal energy stored in the ocean is also unclear, although Maykut and McPhee 
[1995] show that this is important. The spatial distribution of thermal energy in 
the Arctic Ocean on relatively small scales is important in understanding the 
evolution of the sea ice cover. Because water heated by the sun often melts ice, 
which generally allows more heating (the ice-albedo feedback), it is important to 
understand if, where, and when the energy absorbed by the ocean is released. 

The ice-albedo feedback is an example of how spatially varying fluxes in 
the ocean boundary layer can have large-scale, long-time significance. The ice- 
albedo feedback needs to be better understood for accurate modeling of the Arctic 
sea ice cover evolution under climate change scenarios [Curry et al, 1995; 
Perovich et al, 2002a and references therein]. A global climate model (GCM) 
shows the mode of the response to increased greenhouse gas concentrations of 
ocean to ice heat transfer (a heated mixed layer, heated leads, or heated ice) has a 
large influence on ice extent in the Northern Hemisphere [Rind et al, 1995]. 
One-dimensional models also illustrate the importance of ice-albedo feedback. 
Holland et al [1997] show that the ice thickness distribution (particularly the lead 
and first-year ice portion) has a large effect in their model in summer because of 
the solar heating it allows and the resulting ice melt. They report that ignoring ice 
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thickness variability or decoupling the ice thickness distribution from ocean 
conditions will result in an inaccurate response to general model forcing. 

The Arctic influence on the global climate also depends upon the correct 
treatment of the storage and release of thermal energy by the upper ocean [Uttal et 
al, 2002; Randall et al, 1998]. In many global circulation models, increased 
atmospheric CO 2 results in above-average warming in the polar regions due to 
less extensive and thinner sea ice, known as polar amplification (e.g. Manabe et 
al. [1991]). The warming (among other things) varies widely among models due 
to the models’ inability to accurately simulate near-surface temperature structure 
and energy exchange [Randall et al., 1998]. Sea ice feedbacks have been shown 
to play a major part in the global simulations [Rind et al, 1995] and play a part in 
the polar amplification of climate warming trends [Manabe etal, 1991]. 

Changing conditions in the Arctic since about 1990 have added urgency to 
understanding the storage and release of thermal energy in the upper ocean. The 
SEARCH Scientific Steering Committee [2001] describes pan-Arctic change in the 
atmosphere, ice, ocean, and land both over the last decade and over the last 
several decades. Some of the postulated changes most relevant to this study are 
the 43% reduction in mean sea ice thickness [Rothrock et al, 1999], the change in 
atmospheric pressure (mean and spatial pattern) [Walsh et al, 1996; Serreze et al, 
2000], a longer melt season [Smith, 1998], and a warmer, fresher upper ocean 
[McPhee et al, 1998]. Proshutinsky and Johnson [1997] find that two modes of 
wind-driven ice and ocean circulation exist (cyclonic and anticyclonic) in 
numerical studies. The seasonal cycles in these two regimes are very different, 
and Polyakov et al. [1999] argue that the current Arctic climate state is a cyclonic 
regime that in their model causes more summer leads and increased ice export out 
of Fram Strait. The more divergent ice pack allows more open water and 
absorption of solar radiation, potentially giving the ice-albedo feedback more 
influence. The result in the numerical model is a warmer upper ocean, a longer 
melt season, and thinner, weaker, less extensive sea ice cover, in agreement with 
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observations, particularly in the Canadian Basin. In fact, there has even been 
conjecture that an observed shift in Arctic atmospheric circulation could enhance 
melting to the point of a drastically different sea ice regime [McPhee et al, 1998], 

The evolution of the Arctic sea ice cover is not only important for the 
surface energy exchange of the Northern Hemisphere, but also for its effect on the 
world ocean. Sea ice export ultimately freshens the upper waters of the sub- 
Arctic seas. The amount and distribution of fresh water in the sub-Arctic likely 
affects the location and strength of overturning thermohaline circulation in the 
Atlantic Ocean [Aagaard and Carmack, 1989; Carmack, 2000]. The 
thermohaline circulation is a major component of the global heat budget since it is 
responsible for a large fraction of the global poleward heat transport. 

Despite the importance of the ice-albedo feedback, it has been difficult to 
study due to the sparse observations complicated by strong horizontal variability 
in ice conditions and horizontal advection in the ocean. In their study of ocean 
heat storage and release, Maykut and McPhee [1995] were unable to account for 
lateral melting, effects of melt ponds, and lateral advection of water heated far 
away from the experimental site. They were still able to conclude that the 
majority of oceanic heat flux to the ice in the central Arctic comes from solar 
heating of the mixed layer. The Surface Heat Balance of the Arctic Ocean 
(SHEBA) was designed to improve understanding of physical processes 
controlling the vertical and horizontal exchanges of energy in the ocean-ice- 
atmosphere system [Uttal et al, 2002]. The SHEBA experiment was executed 
from the Canadian Coast Guard ship, the Des Groseilliers, which was "docked" in 
an ice floe in the Beaufort Sea in September of 1997. It was allowed to drift with 
the floe until September of 1998. The ship provided a base for an extensive 
program of meteorological, snow and ice, and oceanographic measurements. The 
field program provided an excellent platform on which to carry out a summer lead 
study, both scientifically and logistically. Ultimately, the goal of SHEBA is to 
improve Arctic climate models by improving parameterizations of the physical 
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processes governing surface energy exehange. As illustrated above, without 
detailed observations on small spatial scales, the physical processes sueh as the 
iee-albedo feedbaek would remain a mystery. The present study, while suffering 
some of the previous problems, presents new information and insights on the 
planetary boundary in the summer Arctic. 

The net effeet of the observed melting patterns on the aggregate scale 
(order 100 km) is very difficult to estimate, not only because of strong variability 
of surface conditions on much smaller scales, but beeause of the unknown role of 
mixed layer heating and complex ice-ocean interactions on small scales. It is 
known that the thermal energy of the mixed layer is not strietly dependent on the 
loeal history of lead fraction {Maykut and McPhee, 1995]. The relative motion of 
iee and oeean give rise to a horizontally variable distribution of mixed layer 
temperature often independent of the loeal rate of heating. Since the oceanie heat 
flux to the ice depends both on mixed layer temperature elevation above freezing 
and loeal stress, horizontal gradients in heat flux can be large. Indeed, Wettlaufer 
[1991] observed horizontal variability in Aretie sea iee melt rate (proportional to 
oceanie heat flux) on seales between 10 em and 100 m. Not only is the oceanie 
absorption of solar radiation strongly variable between ice and open water, but 
surface roughness ehanges drastieally among first-year iee, multi-year ice, under- 
iee melt ponds, and leads. 

The small-scale horizontal variability in upper oeean properties and iee 
roughness during the melt season has not been observed until now. Horizontal 
variability has been measured on the top surfaee of the iee, both in regard to 
thermal and mass balance. lee draft has been collected along stake lines and with 
submarine-mounted eeho sounders, but until now no high resolution transeets of 
iee draft were available. Temperature, salinity, and optical properties have been 
eollected in leads, but until now have not been related to the surrounding mixed 
layer. Heat fluxes at discrete depths in the mixed layer have been measured at 
discrete locations on the iee with turbulenee sensors, but the dependence of ocean 
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heat flux on horizontal position and surface type has not been studied until now. 
Vertical profdes of ocean temperature and salinity in the mixed layer have been 
collected during the summer, but tell an incomplete story because they can suffer 
from the unknown effects of advection into different water mass systems and do 
not shed light on mixed layer horizontal structure. All of these improvements in 
studying the ice-ocean system are made possible in this study by use of an 
Autonomous Underwater Vehicle (AUV). 

With new measurement technology come complications and opportunities 
in data interpretation. For the present study, an AUV has been developed for the 
measurement of ice draft, small signals of temperature and salinity and their 
fluxes in the polar mixed layer. The first attempt at taking under-ice horizontal 
profiles of temperature, salinity, and their turbulent fluxes with an AUV was a 
wintertime study by Morison and McPhee [1998]. The concept proved useful, 
and valuable scientific and technical insight were gained. An improved technique 
for estimating turbulent flux from AUV platforms has been developed for this 
study. The new technique and instrument allow collection of a unique data set 
pertaining to the distribution of thermal energy, salt, and even momentum in the 
summertime Arctic mixed layer. 

This dissertation addresses Arctic mixed layer evolution and interaction 
with the surface during the summer melt season with AUV and other SHEBA 
observations and with numerical simulations. Section 2 describes the 
observational study, including methods used to observe the upper ocean and ice 
system, the observations themselves, and a brief discussion of the most important 
results which motivates the numerical studies in section 3. Both one-dimensional 
and two-dimensional models are used to simulate various conditions and to 
compare to observations. Section 4 provides a summary of the SHEBA summer 
story, with input from observations and simulations. Implications of this work on 
larger scales as well as future work and improvements are discussed. 
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2. OBSERVATIONAL STUDY 
2.1 GENERAL DESCRIPTION 

With the success of the Autonomous Conductivity Temperature Vehicle 
(ACTV) in the winter Lead Experiment (March and April 1992; see LeadEx 
Group [1993]) the idea of pursuing the study of spatial variability in the planetary 
boundary layer in the summer case was the logical next step. Very little is known 
of the effect of spatial variability on a boundary layer stabilized by buoyancy flux 
in small patches and bounded by a surface of variable roughness. The solar heat 
flux is also horizontally variable, and the resulting spatial pattern of ice melt is 
unknown. Besides an AUV, other techniques to measure horizontal profiles of 
the under-ice boundary layer do not exist. 

The summer lead sampling strategy is based on previous experience with 
autonomous vehicles in the polar mixed layer. In LeadEx, Morison and McPhee 
[1998] used a complimentary system of fixed and mobile sensors to characterize 
the effects of unstable convection from leads. The set of fixed sensors (turbulence 
instrument clusters or TICs) was placed at a lead’s edge, while the ACTV traveled 
back and forth at several depths under the lead of interest and surrounding ice. 
The ACTV and cluster data showed similar flow statistics when the ACTV was 
upstream and at about the same depth as the cluster. In particular, salinity and 
temperature spectra from the two instruments showed good agreement and a 
wavenumber, k, dependence proportional to The TIC observations were in 
this sense “ground-truth” and allowed the ACTV to accurately map the spatial 
variability of fluxes in areas not upstream of a TIC. 

Figure 2.1 shows the typical summer sampling strategy, along with many 
of the instruments used in this study. Since the buoyancy flux is typically 
stabilizing and therefore suppresses mixing, the depths of the mixed layer and 
fully-developed turbulence are much shallower than in winter. In fact, they are 
often on the order of the deepest ice features, complicating AUV measurements. 
Fortunately, the AUV has an upward-looking altimeter that is useful for mapping 
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ice draft prior to shallow transects. In this study, AUV measurements were 
collected at depths less than 25 m, more often 5-10 m, and for horizontal domains 
approaching 1 km. Once again, fixed sensors were utilized to provide lead-edge 
flow statistics concurrent with cross-lead transects. The spatial variability of ice 
cover was measured with the AUV so that it could be related to the variability of 
the ocean beneath. Numerous Conductivity-Temperature-Depth (CTD) profiles 
near the lead and in the floe interior characterized the vertical ocean structure 
throughout the summer. A Remotely-Operated Vehicle (ROV) collected near¬ 
surface CTD data. Auxiliary measurements like downwelling radiation, ice 
velocity, wind velocity, melt rates, and aerial surveys place the lead 
measurements in context. 

2.2 METHODS 

2.2.1 AUTONOMOUS MICROCONDUCTIVITY TEMPERATURE VEHICLE 
The Autonomous Microconductivity Temperature Vehicle (AMTV) was 
developed to make horizontal profiles of turbulent properties under sea ice. The 
AMTV carries out pre-programmed missions to measure conductivity, 
temperature, depth, ice draft, salinity, turbulent heat and salt flux, magnetic 
heading, and pitch and roll angles and rates. It is based on the REMotely- 
operated Underwater measurement System (REMUS) vehicle developed at the 
Woods Hole Oceanographic Institution. The AMTV is 1.6 m in length, 19 cm in 
diameter, and displaces 33 kg. A photograph is in Hayes and Morison, 2002. 
Pairs of elevator planes and rudder planes at the tail control its motion. The 
AMTV carries a precision Paroscientific pressure sensor and a Systron Donner 
Motion Pack package of accelerometers and sensitive pitch, roll, and yaw rate 
sensors. It employs fast-response Sea-Bird Electronics, SBE 7-02 
microconductivity and SBE 7-01 microtemperature probes, as well as an upward- 
looking Tritech precision acoustic altimeter to measure ice draft. It operates at 
speeds from 1.0 to 1.6 m s'^ under program control. A three-dimensional course 
is programmed into the AMTV. This may be a combination of dead-reckoning. 
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acoustic navigation, and homing to one of two Benthos acoustic transponders. 
The AMTV has the capability to home in on acoustic beacons, either to reach a 
desired checkpoint during a run, or to be recovered in a submerged panel of 
netting. The two-dimensional position relative to a user-defined origin of the 
AMTV is displayed and recorded by a 100-m baseline portable acoustic tracking 
range developed at the Applied Physics Laboratory. In this experiment, the origin 
was the hut location. 

The AMTV has been developed to take advantage of a Kalman smoothing 
scheme that uses vehicle motion to calculate vertical water velocity. The Kalman 
smoother makes an estimate of system state (here, vehicle state) that is an 
optimum combination of a dynamic model-based extrapolation and noisy 
measurements of a few components of the state [Gelb, 1974]. The smoother is 
founded on a linear state space representation of the dynamic system being 
studied. That is, the system evolution is determined by multiple first-order 
differential equations for the rate of change of a vector of state variables. In our 
case, the position and velocity of the AMTV at one time, along with the equations 
of motion allow us to predict its position and velocity at any later time. The 
motion of the vehicle is predicted from the sum of the forces acting on it, 
including buoyancy, the action of water flowing by it vertically, and forces 
exerted by the control surfaces. Because forces on the vehicle and measurements 
have random parts, the system model must include a random forcing function for 
vertical water velocity and an estimate of sensor noise to completely specify the 
system. The smoother uses the vehicle model and noisy observations to produce 
an optimal estimate of the true, unobservable state of the system. In other words, 
the smoother will estimate the portions of the system variability due to unforced 
motion of the body and due to noise. The vehicle forcing (vertical water velocity) 
is calculated to account for what remains. For a complete description of the 
AMTV and the method for estimating vertical water velocity see Hayes and 
Morison, 2002. It is also shown in Hayes and Morison, 2002 how the vertical 
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water velocity spectra from the AMTV and TIC agree over all observed 
wavelengths. 

Covariances of vertical velocity with temperature and salinity are 
calculated to estimate vertical heat and salt fluxes along the vehicle path. A 
“realization” of turbulence quantities is formed by detrending a 100-m segment of 
data (w’, T’, and S’) and calculating the covariance estimates of heat flux, and salt 
flux: 

HF^pCp<WT> (2.1) 

SF = p < w'S'> /lOOO =< w'S'> (2.2) 

The brackets indicate averaging. The 100 m interval is suitable for a realization 
because in that distance the vehicle samples a sufficient number of turbulent 
eddies to produce meaningful flow statistics, yet the distance is not so long that 
the mean conditions change. To determine if the correlation between velocity and 
temperature or salinity is significant, we calculate the correlations at random lags 
and assume they are distributed normally. If the observed correlation falls outside 
of two standard deviations, it is significant [Fleury and Lueck, 1994; Morison and 
McPhee, 199S', Lueck and Walk, 1999]. 

Spectra of vertical water velocity are useful because they allow the 
calculation of two key turbulent parameters: maximum mixing length and friction 
velocity. The AMTV vertical velocity spectra are calculated from 100 m 
segments of data and are averaged in equally spaced bins on a logarithmic 
wavenumber axis. A multi-taper method is used with NW=12 [Percival and 
Walden, 1993]. For each segment, the log of the smoothed weighted spectrum 
[log (kS„„(k)] is fit with a 10* order polynomial as a function of the log of 
wavenumber (log k). 

The maximum mixing length is the vertical scale over which the energy- 
containing turbulent eddies transport momentum, heat and salt. The location of 
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the peak in the spectram of vertical water velocity (in cycles per meter) is taken 
from the 10* order polynomial fit and used to calculate the maximum mixing 
length \McPhee 1994], 


/l = 


0.85 


(2.3) 


The relationship is based on observations that suggest the mixing length does not 
scale with distance from the boundary in the outer boundary layer (outside the 
logarithmic surface layer in the top few meters), but rather with the length scale 
associated with the peak in the wavenumber spectrum of vertical velocity 
\McPhee and Smith, 1976]. 

The friction velocity is a key turbulence parameter representing 
momentum flux 

M, = 4t = \I< u'w'>^ + < v'w'>^ , (2.4) 


where t is the local kinematic stress. Friction velocity is the scale representing 
the overturning speed of the energy-containing eddies. Friction velocity can be 
multiplied by maximum mixing length to produce reasonable estimates of eddy 
viscosity in the entire boundary layer [McPhee and Smith, 1976]. The AMTV 
cannot directly measure friction velocity since the horizontal components of 
velocity are not available. However, an indirect method can be applied as 
follows. 

The spectra of vertical water velocities collected by the AMTV can be 
used in a novel calculation of friction velocity using a technique developed by 
McPhee [2003], which is recapitulated in some detail in Appendix A. If the 
spectrum contains a wavenumber range where the turbulent kinetic energy 
depends only on production of energy and wavenumber of disturbance (an inertial 



13 


subrange), then it can be shown that there exists a universal, non-dimensional 
spectrum in this range. The friction velocity is calculated as the value that 
nondimensionalizes the observed spectrum to the universal spectrum. This 
technique is similar to the inertial dissipation method {Edson et al, 1991] in that it 
uses spectral densities of turbulent vertical velocity to calculate momentum flux 
(stress). However it does not require the simplified turbulent kinetic energy 
equation of production equals dissipation, nor does it require specification of the 
Kolmogorov constant. Rather, simple dimensional analysis in the inertial 
subrange is used. The calculation is described in Appendix A. 

2.2.2 SETTING AND AUXILIARY MEASUREMENTS 

The present study relies heavily on data gathered during a SHEBA 
intensive summer observation period from July 10 to August 10, 1998. The ship 
track and ocean bathymetry for the SHEBA drift period are shown in Fig. 2.2. 
June, July, and August are the blue, pink, and black segments near the end of the 
drift. At the beginning of the summer, the ship drifted from the Beaufort Sea onto 
the Chukchi Cap, where it drifted slowly until late July, when a rapid drift off the 
Cap began. 

During the intensive summer observation period a suite of measurements 
was taken near “Sarah's Lake”, a lead that opened in late-May in an area of first 
year ice previously used as the airstrip. The base of operations was a hut located 
about 1 km off the port bow of the Des Groseilliers, near the lead edge. Outlines 
of the lead at different times are shown in Fig. 2.3. The perimeter of the lead was 
traversed with a small boat towing an acoustic beacon that was being tracked. 
The lead remained approximately the same size and shape for most of July, and it 
closed in the first few days of August. It opened to a much larger size by Aug. 4 
(day 217), and even larger size by day 219 (not shown). 

The basic oceanographic measurements took place at "ocean city" which 
was in nearly continuous operation throughout the SHEBA year. A continuously- 
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operating CTD (SBE 911+) collected profiles to 150 m several times per day with 
the exception of a few instrumental problems in July. 

Data on ice-ocean exchange of momentum, heat, and salt were collected at 
ocean city, and they provide a valuable comparison with the lead measurements. 
From mid-July to mid-August, the ocean city TIC was operating between days 
211 and 215 and for days 219-221. Low ice speeds and biofouling led to the data 
dropouts. Turbulent Instrument Clusters (TICs) consist of a rigid stainless steel 
mast holding two instrument clusters at set depths (both less than 10 m in our 
case). Each cluster consists of a set of orthogonal partially ducted current meters, 
a standard SBE-04 conductivity cell or SBE 07-2 microconductivity probe, a SBE 
07-1 microtemperature probe, and a pressure sensor to collect mean quantities of 
velocity, temperature, salinity, and depth. The microtemperature and 
microconductivity probes have a fast response, sampling at 6 Hz and averaging 
over a 1- to 4-second period prior to recording to computer disk. The instruments 
allow the calculation of turbulent quantities of variance and covariance in the 
velocity components, temperature and salinity. In particular, we are interested in 
the vertical fluxes of momentum, heat, and salt, which are calculated using Eqs. 
2.1, 2.2, and 2.4. An averaging interval of 15 minutes is used in forming the 
covariances, analogous to the 100 m AMTV segments. The speed of the ice 
relative to the ocean at the cluster depths must be greater than about 5 cm s'' to 
overcome the rotor threshold. McPhee [1989] gives a more complete description 
of the TIC, and McPhee [1992] establishes the TIC as a reliable system for 
measuring turbulent quantities in the under-ice boundary layer. 

At the lead site, a TIC using SBE-03 temperature probes and SBE-04 
conductivity probes collected data near the lead edge. A cluster was typically 
located at the same depth as the AMTV operating depth. Flow statistics from the 
TIC could be compared to upstream AMTV segments. The relative ice-ocean 
velocity at the cluster depths did not generally exceed the threshold until day 205, 
corresponding to a small peak in wind velocity, then again on 208 when the wind 
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began to increase. The TIC was damaged by a ridging event on day 208. On 
days 205 and 208 there was a cluster at 3.5 m and 5.5 m from the surface. The 
lead TIC was resurrected on days 219-220, with a cluster at about 5.5 m. In the 
latter period the water traversed most of the lead before reaching the TIC. 

Also at the lead site, we operated a SBE-19 CTD in mooring and profding 
modes to help describe the local temporal and spatial variability during the 
summer observation period. The vertical profiles were useful in field calibrations 
of the AMTV temperature and conductivity probes, as well as numerical model 
initialization. They provided profiles more representative of conditions in the 
lead than the ocean city profiles, which were collected almost one kilometer 
away. The CTD was suspended from the ice at AMTV operating depths to serve 
as a check on the advective variability. 

Yet another SBE-19 CTD was attached to a small Remotely Operated 
Vehicle (ROV) to sample the upper few meters of the ocean under various surface 
types. The ROV had a video link to the hut, along with a real-time display of the 
CTD data. Visual ice properties, ice draft and surrounding water properties were 
observed. These data were saved to videotape and computer disk. We obtained 
many short vertical profiles with the ROV under the ice and lead by driving it to 
10-12 m, then allowing it to float to the surface. The ROV was operated from the 
same hut as the AMTV. The ROV was also used to aid in or observe the recovery 
of the AMTV after each mission. 

Auxiliary measurements included ice velocity, wind velocity, 
downwelling solar radiation, area coverage by surface type (ice, melt pond, lead), 
lead albedo, ice and melt pond transmittances, surface, lateral, and basal melt 
rates, ice permeability, large-scale ice deformation, and a host of others related to 
atmospheric processes. 

Ice velocities are based on Global Positioning System (GPS) 
measurements made at the ocean city site. In the calculations and plots here, 
inertial motion has been removed with a complex demodulation algorithm 
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\McPhee, 1986a, Morison et al, 1987, McPhee et al, 1987]. It has been shown 
that the interfacial stress due to inertial oscillations is generally not significant, 
i.e., the ice and mixed layer tend to oscillate together [McPhee, 1986b]. Mean 
velocity is calculated at a time interval of three hours. 

Ten-meter wind velocity and downwelling broadband “shortwave” 
radiation were measured by the SHEBA Project Office meteorological tower and 
averaged over one hour. The shortwave detector was an Eppley PSP with a range 
of 285-2800 nm. The amount of surface area covered by ice, leads, and melt 
ponds was measured from photographs collected on helicopter surveys throughout 
the summer [Perovich et al, 2002b]. The surveys were roughly every week from 
May to October of 1998, and the flight path was a 50 km by 50 km box centered 
on the ship. Full coverage flights on a 10 km by 10 km box centered on the ship 
were conducted on 17 May and 25 July. A few radiation measurements were 
made under ice for the purpose of spectral transmittance calculations [Light et al, 
2003]. The albedo of the lead at the SHEBA site was measured [Pegau and 
Paulson, 2001]. The changes in thickness and extent of the ice pack due to 
surface, basal, and lateral melting were observed at several sites and are 
summarized by Perovich et al [2003]. Many types of ice (ridges, ponds, first- 
year, multi-year, etc.) were instrumented with ablation stakes and thickness 
gauges and thermistor strings. 

2.3 OBSERVATIONAL RESULTS 

The early and middle summer periods at SHEBA (mid-June to late-July 
1998) were characterized by low wind and ice velocity, strong downward 
radiative fluxes, and intense surface melt that collected in leads and ponds. The 
snow depth averaged only 3 cm by the end of June. In late July and early August, 
high wind and ice velocities rapidly mixed the fresh surface waters downward and 
rearranged the ice cover. Surface melt tapered off, and basal melt increased along 
with lead fraction. By late August, the leads began to freeze. 
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2.3.1 GENERAL RESULTS 

Wind stress and melting control the behavior of the ice pack over the 
su mm er season. The atmosphere-ocean interface in summer leads is 
characterized by strong solar radiation, less stress than the ice-ocean interface, 
and pooling of fresh water due to lateral flow of meltwater. Here we examine the 
observed exchanges of thermal energy, fresh water, and momentum across the 
various types of upper boundaries of the ocean. We begin with the bulk exchange 
of properties. 

The surface friction velocity, u*o, is a defining characteristic of the 
planetary boundary layer and is shown for SHEBA in Fig. 2.4a. The surface 
friction velocity is defined in terms of the surface stress, 

Tq = , (2.5) 

where p is the water density. For Fig. 2.4a, the surface stress is estimated from 
the 10-m wind speed, Uio, and ice speed Uu according to Eq. 2.6 \McPhee, 1990]. 
(Bold indicates a vector quantity.) This assumes that internal stress gradients in 
the ice are negligible, a good assumption during summer in the Arctic when the 
percentage of open water is increased. Steady state ice motion is assumed; an 
assumption that may be poor when the wind changes quickly over an inertial 
period. Horizontal homogeneity is also assumed (open water is not accounted 
for). 

To = PuCdo^ioUxo ( 2 - 6 ) 

The density of air, pa, is 1.3 kg m'^; the drag coefficient, Coa, is 1.5x10'^. The 
drag coefficient was shown to vary from 1.1 to 2.0x10' during SHEBA [Andreas 
et al, 2001]. This range results in upper and lower bounds on friction velocities 
during the storm events of ±15% (less during quiet periods). Density (/?,) and 
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thickness {hi) of ice are approximately 920 kg m' and 1 m. The Coriolis 
parameter,/ is approximately l.dxlO'^ s'' at 80° N. 

The heat available in the mixed layer to melt ice is proportional to 8T, the 
elevation of the temperature, T, above the surface freezing temperature, T/. 

Sr = T-Tf (2.7) 

The elevation of the observed 4-m temperature above freezing is shown in Fig. 
2.4b, where Tj- » -mS and m = 0.055 "C psu~^. 

The rate at which heat energy in the mixed layer reaches the ice, the heat 
flux, can be related to the temperature elevation and the surface stress through an 
empirical bulk relation. 

<wT')o = CfjpCpU^^ST (2.8) 

The exchange coefficient, ch is approximately 0.006 over a wide range of 
conditions [McPhee, 1992]; Cp is the heat capacity of seawater at constant 
pressure (3980 J kg''). Figure 2.4c shows the ice-ocean heat flux (positive 
upward), calculated from the bulk exchange formula. Heat flux was upward at all 
times because the mixed layer was always above the freezing temperature. The 
solid line shows heat flux using a mixed layer temperature that has been 
interpolated to fill in missing data. The circles show heat flux when temperature 
elevation was actually observed. 

Figure 2.4d shows the bulk estimate of salt flux at the ice-ocean interface 
that results from the ocean heat flux of Fig. 2.4c. Upward heat flux at the ice- 
ocean interface, if not compensated by conduction through the ice, causes ice 
melt. The rate of release of meltwater downward into the water column (upward 
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salt flux) due to the oceanie heat flux can be estimated using the salt balance 
equation \McPhee, 1992], 

{WS'),={w'r),{S„,-S,Jq~^ (2.9) 

For Fig. 2.4d, the ice salinity is taken to be 2 psu (varies from 0 to 4 psu 
depending on ice type and vertical position according to Sicken et al. [2002]), and 
the value of qi is set at a representative value of 74 K. The assumption that there 
is no heat conduction from the ice to the atmosphere is reasonable because most 
of the ice is isothermal at the freezing temperature for most of the summer 
\Perovich et al, 2003]. However, it is important to note that meltwater can also 
enter the ocean via runoff or percolation from surface melting, but it is not 
included here. 

Ocean city CTD data for June, July, and August are presented in Fig. 2.5. 
Bathymetry during the SHEBA drift is shown with salinity contoured in depth and 
time (Fig. 2.5a). In Fig. 2.5b, the square of the buoyancy frequency, N , is 
contoured. 

( 2 . 10 ) 

p 8z 

where g is the gravitational constant. The temperature is not used in the 
calculation of density gradient here because it affects the density very little near 
the freezing point. An objective, if arbitrary, definition of the mixed layer is 
arrived at by finding the depth at which the squared buoyancy frequency exceeds 
a specified value. The 1x10' s' contour is plotted in magenta and taken to be the 
mixed layer depth. (Choosing a value of 5x10'"^ s'^ would lower the mixed layer 1 
to 2 m.) Also shown in Fig. 2.5b are times and depths of AMTV runs (black 
circles). After day 210, most of the runs were in the mixed layer. Potential 
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temperature and elevation of potential temperature above surfaee freezing 
temperature are also shown (Figs. 2.5c and 2.5d). 

2.3.2 SEQUENCE OF EVENTS WITH SPECIFIC AMTV EXAMPLES 

The SHEBA summer experiment can be divided into two periods 
separated by a transitional period. The first period of mid-June to late-July was 
relatively calm: Fig. 2.4a shows the surface stress and wind were low until a 
moderate peak centered on day 208, followed by a large, sustained event from 
about 210-214 (the transitional period). The ice-ocean heat flux was also 
relatively weak before day 208 (Fig. 2.4c). The vertical structure of salinity from 
ocean city (Fig. 2.5a) changed little during the quiet period, with the exception of 
the period when the camp drifted onto the Chukchi Cap (days 159-165). In 
addition, the mixed layer shoaled from about 15 m to 6 m (Fig. 2.5b). 

The mixed layer warmed in the pre-storm period, as shown by potential 
temperature (Fig. 2.5c) and elevation of potential temperature above surface 
freezing temperature (Fig. 2.5d). Other temperature features during the quiet 
period were concurrent with the drift onto Chukchi Cap. 

In the pre-storm period, snow and ice melt was trapped near the surface in 
leads, cracks, and ponds. The small-scale vertical structure of salinity and 
temperature near the lead was measured with the ROV CTD and a CTD towed 
from a small boat [Paulson and Pegau, 2001]. Paulson and Pegau [2001] found 
that a warm, fresh layer of meltwater grew steadily in the SHEBA lead from mid- 
June to late-July and shrank rapidly during the main storm on days 210-214. The 
very fresh surface layer is not visible in the ocean city profiler data in Fig. 2.5 
because the ice around the hole was thicker than the fresh layer. The ROV 
camera and CTD tracked the depth progression of the fresh-water interface in the 
SHEBA lead (Fig. 2.6a and 2.6b). On day 206 (July 25) the layer reached a 
maximum depth of 1.2 m in the lead while the ice draft surrounding the lead was 
slightly less. The layer had a salinity of about 2 psu and a temperature of 1.6 °C 
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[Paulson and Pegau, 2001]. The ROV was unable to sample the uppermost part 
of the layer beeause it was insufficiently buoyant in the low-density water there. 
Leads sampled within 30 km of SHEBA on day 203 showed similar warm, fresh 
layers around 1 m thick [Richter-Menge et ah, 2001], similar to observations from 
the SHEBA lead. The fresh water formed from surface melt did not appear to mix 
into the upper ocean during the quiet period. 

The AMTV data collected during the quiet period (day 204) indicate that 
vertical mixing of temperature and momentum was weak. Figure 2.7 shows 
vehicle depth and ice draft in the first panel, temperature in the second and 
horizontal position in the third. The vehicle was at a depth of 10 m initially; it 
climbed to 6 m then back to 10 m (segments 1-3, respectively). The mean ice 
draft was just over 1 m, with a ridge of about 3 m draft. The temperature profile 
indicates a weak vertical temperature gradient: temperatures at 6 m (segment 2) 
were about 0.02 C° lower than at 10 m. The outline of the lead is shown in the 
third panel. The star represents the starting point of the run, the circles and 
squares indicate the beginning and end of each run segment to be analyzed. On 
these “legs” the AMTV was not making maneuvers such as depth or heading 
changes, allowing one to obtain vertical water velocity using the Kalman 
smoother. The ice velocity of 6.4 cm s ' (from demodulated GPS data) is shown 
as an arrow; the length of the arrow corresponds to the distance the ice traveled in 
1000 sec (about 17 min). The thick line segments indicate open water acoustic 
returns. This run is typical of runs made from day 203 to day 208, as the lead 
geometry and direction of ice motion were fairly constant. 

Figure 2.8 shows the deviation quantities of detrended temperature, 
vertical water velocity, heat flux, and friction speed. Unfortunately, the 
microconductivity sensor was damaged slightly in a recovery operation on day 
195 and was not operating correctly. It was replaced on day 213. No salinity- 
related quantities at the lead are available from the AMTV from day 195 to 212. 
Temperature shows a smaller variance at 6 m (segment 2) than at 10 m (segments 
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1 and 3), and no effects of the lead are visible. Vertical water velocity 
fluctuations are fairly uniform with depth and horizontal position, with perhaps 
more variance at 10 m. The covariances of temperature and vertical water 
velocity are calculated from 100-m segments of data. A horizontal bar represents 
the value of each segment. Only the first two bins downstream of the lead (the 
300 m to 500 m region in the distance coordinate) at 6 m depth are statistically 
significant at the 95% level. The fluxes there are 10 W m'^ upward followed by 3 
W m'^ downward. The friction speed is shown with the surface friction speed of 
Fig. 2.4. The small dots represent the 95% confidence interval. Friction speed is 
steady at 0.5 cm s'' except for one point near the end of the run. Mixing length is 
steady at about 1.2 m throughout the run (not shown). There were relatively low 
friction speeds and negligible or small heat fluxes at 6 and 10 m near the lead on 
day 204 (during the quiet period). 

Large changes were observed in the ocean during the first major storm of 
the summer around days 210-214. The wind reached its highest velocity of the 
previous few weeks (Fig. 2.4a). The ice-ocean relative velocity increased enough 
around day 210 to detach long diatom strands from the ice underside. The strands 
had been growing under the ice near the lead and had reached lengths of up to 1 
m. The release of the mucus-like strands was problematic for many instruments. 
(From day 206.5 to day 208.5 as the wind began to strengthen, the ocean city 
CTD was not functioning primarily due to biofouling.) The increase in wind 
speed was associated with an increase in bulk ice-ocean heat flux. The heat flux 
was the strongest of the summer from day 208 to day 212 (Fig. 2.4c). Although 
the wind event around day 208 was smaller than what followed, the mixed layer 
contained more thermal energy, which gave rise to a comparable flux magnitude. 
By the end of the storm, the mixed layer cooled to about 0.1 C° above the freezing 
temperature (Figs. 2.4b and 2.5d) and deepened to about 16 m (Fig. 2.5b). A 
layer of warm water between 20 and 30 m in the pycnocline remained through the 
storm (Fig. 2.5c). Most dramatically, from day 209-214, the mixed layer 
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freshened more than 1 psu (Fig. 2.5a). As the camp drifted from Chukchi Cap 
into deep water during the storm, the pycnocline also freshened from 30-60 m 
(Fig. 2.5a). 

The change in mixed layer salinity during the storm can be used to 
estimate the total salt flux at the ocean interface. The ocean city salinity profdes 
in Fig. 2.5a are used to calculate the rate of change in salt content of the top 30 m. 
This may include some advective changes not due to surface processes, especially 
during periods when the camp drifted over topographic features. The advective 
rate of change in integrated salt content is that estimated by assuming the whole 
water column salinity changed by an amount represented by the salinity change at 
30 m. This is subtracted from the actual rate of change of salt content in the top 
30 m to obtain the interfacial flux. We find a significant flux only from days 210 
to 215: 4.6 ± 1.5x10' kg m' s', corresponding to about 80 cm of fresh water 
mixing down over those 6 days. The 95% confidence interval is formed from the 
uncertainties in the linear best-fit model used to calculate the rates of change. We 
find a similar result with the CTD data from the lead edge. Because the salt flux 
includes meltwater stored in leads and percolating through the ice, it is about an 
order of magnitude larger than the instantaneous salt flux due to bottom melt 
shown in Fig. 2.4d. 

The transition period is largely signified by this mixing away of the fresh 
layer of water in the lead. During the first wind event from 206.5-210.5, the layer 
receded from 120 cm to 65 cm (Fig. 2.6a and 2.6b). This is equivalent to a salt 
flux in the lead of 5.0± 1.0x1 O' kg m' s' . During the major wind event from day 
210.5 to day 212.5, the remaining 65 cm of fresh water were mixed away, 
corresponding to a salt flux of about 12x10'^ kg m'^ s'\ The 50 km x 50 km area 
around SHEBA began to diverge around day 209 [Stern and Moritz, 2002], but 
the SHEBA lead itself closed by day 213 before widening in the days following. 
The effect of divergence or convergence on the thickness of the fresh water layer 
is neither known nor accounted for in the salt flux calculation, however it is 
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possible that divergence could thin the fresh layer. In this case a vertical salt flux 
out of the lead is not implied. 

The ROV video record clearly revealed the change under the ice as the 
fresh layer eroded. The diatom colonies were ripped from the ice underside, and 
ice crystals formed in the top few meters of the ocean as the fresh surface water 
interacted with the colder, saltier water on day 209. The frazil ice floated to the 
surface and encased what remained of the diatom colonies, before eventually 
melting away a day or so later. 

An AMTV run from the storm period shows increased activity in the 
mixed layer: large momentum, heat, and salt fluxes were observed. Figure 2.9 
illustrates the run setup. On day 213 (August 1), the major wind event of the 
su mm er was in full swing, with wind speeds of over 8 m s'^ (16 knots). The ice 
cover shifted around changing the geometry of the lead. Open water sections 
were no wider than about 50 m (Fig. 2.9a). In Fig. 2.9c, we see the AMTV 
initially traveled in the negative x-direction (segment 1) and returned past the 
origin to the positive x direction (segment 2), returning (segment 3) to make a 
final excursion in the negative y direction (segments 4 and 5). The ice velocity is 
shown: 21 cm s '. The AMTV was at a depth of 8 m, and observed a mean ice 
draft of about 1 m with a ridge approaching 4 m depth located at x=150 m, y=0 m 
(sampled in the last half of segment 2 and halfway through segment 3). 

Deviatory data (now including salinity deviations) are plotted in Fig. 2.10. 
Absolute salinity is questionable due to a nonlinear drift in the microconductivity, 
and the first 500 m of conductivity data were unusable. The AMTV encountered 
large temperature variations as it traveled from -250 m to +250 m in the x- 
coordinate and returned to the origin (segments 2 and 3 in Fig. 2.10a— 
particularly near the ridge). Plots of salinity deviations and vertical water velocity 
show only slightly elevated activity in this region (Fig. 2.10b and 2.10c). Note 
the vertical water velocity amplitude was roughly twice that of day 204. Figure 
2.11 shows significant turbulent fluxes from 500 m to 1100 m in the run 



25 


coordinate (the portions of segment 2 and 3 under and downstream of the ridge as 
indicated). The heat flux was 72-74 W m' upward for three 100-m segments and 
31 W m'^ upward for one. Salt flux also had significance for four segments: one 
was 0.82 kg m'^ s'* upward and three of them averaged 1.1 kg m'^ s'* downward, 
an apparent paradox with the expected upward heat flux and upward salt flux 
under melting ice. Friction speed was high, especially for segments under and 
downstream of the ridge. Values were even above the surface value of about 1 
cm s'*. The friction speed was 2-3 times larger than on day 204. Maximum 
mixing length from vertical water velocity spectra was steady at 1.0 m, except for 
two extrema at 1.5 m and one at 0.2 m. The AMTV observed an active mixed 
layer on day 213 during the storm, particularly under and downstream of a ridge. 

After the storm a brief quiet period around day 214-217 was followed by 
another wind event from 218-220 (Fig. 2.4a). Combined with a brief rise in 
temperature, this wind event corresponded to a smaller, briefer pulse of heat flux 
around day 218 (Fig. 2.4c). The mixed layer salinity and temperature were 
otherwise relatively steady (Figs. 2.5a and 2.5c), and the mixed layer depth was 
around 20 m except for a brief shoaling to 10 m (Fig. 2.5b). The layer of warm 
water in the upper pycnocline was maintained following the major storm and 
throughout the rest of the summer: the mixed layer was within 0.1 C° of the 
freezing temperature, while the pycnocline was about 0.2 C° above freezing (Fig. 
2.5d). Another mid-depth temperature maximum was observed over the sloping 
bathymetry (Fig. 2.5c) as camp continued to drift off Chukchi Cap, presumably 
related to topographically steered currents along the shelf slope. 

In the post-storm period, meltwater was not trapped near the surface but 
entered the ocean as the melting occurred. The fresh layer of meltwater did not 
reappear. The lead profiles of temperature and salinity generally resemble the day 
219 (August 7) profile in Fig. 2.6. A detailed view (Fig. 2.6c and 2.6d) shows the 
lead surface had weak gradients, was stably stratified, and was above the freezing 
temperature. 
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For the post-storm period, the AMTV data indicate strong fluxes of 
momentum in the lead and under ice, but only strong heat and salt fluxes in the 
lead. The last day of AMTV observations was day 219-220 UT. The wind had 
increased again after a few days of relative calm, and the lead had grown 
dramatically in size around day 215. Figure 2.12 shows the run on day 219.9. 
Vehicle depth was 5 m, and ice draft averaged about 1 m with several ridges of 
about 3 m (Fig. 2.12a). After traveling for 800 m without leaving the lead 
(segment 1), the AMTV returned (segment 2) and continued to travel in the 
positive x-direction 200 m under the ice (segment 3) before returning to the origin 
(segment 4). See Fig. 2.12d. The ice velocity of 15 cm s ' is shown. The 
temperature increased as the lead was approached from upstream (segment 2), 
peaked at the lead edge, then fell with distance under the ice (segment 3). See 
Fig. 2.12b. Nearly the same trends were observed on both the “out” and “back” 
legs. Salinity decreased over segment 2 then remained uniform under the ice. 
The overall slope of salinity (particularly segment 1) is biased due to the 
nonlinear drift in the conductivity. Figure 2.13 shows both temperature and 
salinity had much higher variances in the lead compared to under the ice. Also, 
temperature and salinity deviations were about four times bigger than on day 213, 
while vertical water velocity was slightly less energetic than day 213. Figure 
2.14a and 2.14b show strong downward heat fluxes and upward salt fluxes under 
the lead and smaller fluxes under the ice. The lead fluxes of heat and salt 
averaged 98 and 144 W m'^ and 1.7 and 2.4x10'^ kg m'^ s"' (segments one and 
two). Segment three did not pass the significance test for heat or salt flux, and 
segment four only had significant salt flux: 1.9x10'^ kg m'^ s’'. The segment 
averages reported here are calculated from the covariances over each entire 
segment, not the average of the 100 m segments. The first 200 m of segment one 
were not included due to salinity drift contamination (see Fig. 2.13). Friction 
speed was close to the surface value and between the friction speeds of days 204 
and 213 (Fig. 2.14c). Maximum mixing length averaged around 0.7 m with a 
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maximum of 1.8 m at the lead edge. AMTV data from after the major storm of 
the SHEBA summer show exchange at the ocean surface, especially in leads. 

2.3.3 BULK, TIC, AND AMTV OBSERVATIONS 

The bulk estimates of fluxes of Figure 2.4 are compared to the values just 
below the ice measured directly by the ocean city TIC in Fig. 2.15. The TIC 
measurement depths were 3.5-4 m (cluster 1) and 7.5-8 m (cluster 2), and the ice 
draft at ocean city averaged 1.4 m from day 210 to 220. Figure 2.15 shows the 
friction speed, heat flux, and salt flux from clusters 1 and 2 for the period along 
with the bulk estimates. Each data point represents an average of several 15- 
minute “realizations” of turbulent quantities (from 4 to 12 realizations per data 
point). The salt flux estimates from the deeper cluster were corrected for the 
temporal lag of the standard SBE-04 conductivity sensor [McPhee and Stanton, 
1996] with output of a microconductivity cell. 

Fluxes below the ice can be compared directly to the surface estimates and 
to a similarity prediction based on the surface values. The solid points in Fig. 
2.15 are predicted for the TIC depths from the analytical similarity solution for 
the horizontally homogeneous, steady boundary layer forced by surface stress and 
buoyancy flux [McPhee, 1983]. The surface stress is determined from force 
balance, and the buoyancy flux (proportional to salt flux) only includes melt 
estimated from ice-ocean heat flux—^no runoff. The theory predicts an 
exponential fall off with non-dimensional depth, of momentum and scalar 
fluxes. 

4 ,. ( 2 . 11 ) 

( 2 . 12 ) 


(2.13) 
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where H(f) = e*«, = ^ = z(^)-', r,:^ ^i+ M"'Oo 


w'b')^=gas{wS')^ -ga^{w'T 


(2.14) 

(2.15) 


The last equation is the buoyancy flux, and as and aj-are expansion coefficients 
for salinity and temperature. Constants are a:= 0.4 (the von Karman constant), 
4)v=0.05, A*=0.028, and Rc=0.2 (the critical flux Richardson number). See below 
for more discussion on length scales in turbulent boundary layers. The surface 
buoyancy flux enters the calculation through the stability parameter, rj*. 

The friction speed and fluxes observed by the TIC are significantly lower 
than the surface values and those predicted by similarity theory for days 211-215. 
On days 219-220 the TIC friction speeds are close to the surface value and even 
higher at times (Fig. 2.15a). The lower and upper clusters measured similar 
stresses on all days. TIC heat fluxes are smaller than expected for days 219-220, 
and the lower cluster measured weaker fluxes (Fig. 2.15b). Contrary to the earlier 
period, salt fluxes on day 219-220 are much higher than estimated at the surface 
from the heat balance (Fig. 2.15c). Recall that the salt flux at the upper ocean 
boundary based on mixed layer salinity changes from day 210 to day 215 is about 
10 times larger than the heat balance estimate shown in the solid blue line of Fig. 
2.15c. Using the larger flux gives more strongly attenuated similarity predictions 
for stress and heat flux, which is addressed in numerical simulations to follow. 

The AMTV run averages compare favorably to the surface values of Fig. 
2.4 and lead-edge TIC values. The friction speed comparison is shown in Fig. 
2.16. The blue line indicates values calculated from wind and ice velocities as in 
Fig. 2.4. While the AMTV-derived friction speeds vary greatly over any given 
run, the run averages of Fig. 2.16 agree reasonably well with the surface estimate. 
The average is computed from the 100-m segment values like those of Fig. 2.14. 
Pink crosses show the averages for each run. On many days, especially before 
day 209, the AMTV sampled multiple depths in the same run. Segments within a 
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run that were above a non-dimensional depth of 0.4 (see Eq. 2.14) based on bulk 
surface estimates (no percolation) were averaged to form the red squares. With 
only two exceptions this amounts to averaging segments that were above 10 m 
(see Fig. 2.5b for run depths). These points are considered to be in the boundary 
layer. Agreement with the surface friction speed is good, except for more 
extreme values of the peaks and valley on days 211-213. Lead-edge TIC data is 
shown with the blue stars (3.5 m depth) and green stars (5.5 m depth). Except for 
day 205, the TIC shows friction speeds at or above the surface value. On days 
219-220, the lead edge TIC (only one at 5.5 m) mostly agrees with the AMTV. 
This can be contrasted to the ocean city TIC (Fig. 2.15) where the observed 
friction speed was typically less than the surface estimate except on day 219-220. 
On day 219-220 both TIC systems and the AMTV were operating, and there is 
considerable overlap. 

Figure 2.17 is a similar AMTV plot for heat flux. The heat flux calculated 
from the bulk exchange formula is the solid blue line (as in Fig. 2.15). AMTV 
run averages for all depths and for depths in the boundary layer are shown by the 
pink crosses and red squares, as in Fig. 2.16. Only those 100 m segments that 
show a significant correlation between vertical water velocity and temperature are 
used in the averages. The AMTV averages are mostly upward but not as strong as 
the bulk estimate. One exception is day 213, where heat flux is much larger 
upward. Other exceptions are the data after days 215 when the AMTV heat flux 
is strong and downward. The ice fraction increased dramatically after day 215, so 
most of these segments are not under ice. The TIC heat flux from the downstream 
edge of the lead is always downward at both 3.5 m and 5.5 m depths. The ocean 
city TIC heat flux (Fig. 2.15b) is not comparable to the lead-edge TIC, but does 
show a suggestion of the high upward heat flux on day 213, followed by lower 
values on day 215, like the AMTV. However, at no time does the ocean city TIC 
indicate downward heat flux. 
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We suspect the strong downward heat flux at the lead after day 215 is due 
to radiative heating. Figure 2.18 shows the detail on days 219-220, when the 
downward heat flux at 5 m in the lead was especially strong. The additional solid 
blue line is the downward shortwave radiative flux as measured by the SHEBA 
Project Office. The lead TIC and AMTV are in fair agreement, and the TIC heat 
flux at 5.5 m is consistent with the shortwave flux, at 30-50% of the amplitude. 

The AMTV average salt flux comparison is illustrated in Fig. 2.19. As 
mentioned above, the AMTV conductivity probe was not functional until day 213. 
Until day 219, very few AMTV runs showed significant correlation between 
vertical water velocity and salinity. The lead TIC showed upward salt fluxes on 
days 205 and 208, and those were especially large on day 208 when the first wind 
event was peaking. The AMTV salt flux on day 213 was downward, as described 
(Fig. 2.11). A strong upward salt flux was observed on day 215, and then small 
values prevailed until 219. For days 219-220, the AMTV and lead TIC both 
showed strong upward salt flux well above the surface heat balance estimate, as 
did the ocean city TIC on day 219 (Fig. 2.15c). 

2.4 DISCUSSION 

2.4.1 PERIOD I: SURFACE TRAPS MELTWATER 

For much of the melt season, ice and snow meltwater was trapped near the 
surface. From day 152-190 (June to early-July), the mixed layer temperature and 
salinity were relatively steady (Fig. 2.5a and 2.5c), yet the surface was rapidly 
melting \Perovich et al, 2003]. Since the ice was still not permeable, surface 
melt drained laterally into low-lying areas of melt ponds and leads [Eicken et al, 
2002]. Mid-July (days 190-207) saw that the mixed layer warmed in the top 30 m 
and shoaled but still did not freshen significantly. Higher ice permeability 
allowed more vertical percolation of melt in mid-July, but fresh water continued 
to build up in leads and cracks [Eicken et al, 2002; Paulson and Pegau, 2001]. 
The buoyant melt probably tended to find the sea level (the lead surface) even 
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when it percolated through the ice underside. This means lateral flow of 
meltwater was still significant when vertical percolation was occurring. 

There are two reasons why surface melt did not mix into the upper ocean 
during the quiet period. First of all, fresh water that entered the ocean often 
formed an underwater ice layer or false bottom when it came into contact with the 
colder seawater: false bottoms tended to trap the runoff near the surface [Eicken 
et al, 2002]. Underwater ice forms when fresh water and seawater at or near their 
respective salinity-determined freezing temperatures contact each other under 
quiescent conditions {Untersteiner, 1961; Martin and Kauffman, 1974]. Because 
solar heating was larger in leads, the meltwater that was collected there and the 
seawater beneath it were above their respective freezing temperatures, and 
underwater ice could not form. The second reason fresh water did not mix 
vertically was because of the strong and stable density gradient near the lead 
surface and low ice speeds (typically 5 to 7 cm s'' in mid-July). As described 
above and in Paulson and Pegau [2001], meltwater intermittently formed a fresh 
surface layer in the SHEBA lead from day 170 to day 184 and grew steadily to 
1.2 m by day 207. 

It is possible to quantify the tendency for a buoyant layer in a lead to form 
with an Obukhov length, Lq. 


L 


0 


K{w'b\ 


(2.16) 


The Obukhov length is the distance from a boundary at which the buoyancy 
production or destruction of turbulent kinetic energy is balanced by shear 
production. At depths less than \Lo\, turbulence is generated by shear, and 
buoyancy effects are small, which is known as forced convection. In wintertime 
unstable or free convection, dense brine rejected during ice formation drives 
cellular convection and generates turbulence for depths greater than \Lo\. In the 
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summer, stabilizing meltwater tends to balance shear production. During periods 
of freshwater flux, mixing (shear-driven) is confined to depths shallower than \Lo\. 
During times of relatively high wind stress, the mixed layer may be shallower 
than \Lo\, and the mixed layer would tend to deepen. 

An estimate of the Obukhov length beneath the lead can be made for days 
190-207. The surface friction speed rarely exceeded 0.005 m s ' and averaged 
0.0035 m s”' (Fig. 2.4a). An upper bound on the freshwater flux in the lead is 
found from the mean rate of surface ablation from day 190-207: 1.5 cm day'* 
[Perovich et al, 2003]. Assuming all of the melt from ice drained into leads 
(about 75% and 5% of area, respectively, from Perovich et al. [2002b]) gives an 
equivalent melt rate of about 23 cm day *. This rate of ice melt corresponds to a 
freshwater runoff flux into the lead by Eq. 2.17 of about 7x10'^ kg m'^ s'* or 
5.3x10'^ W kg'* for buoyancy flux by Eq. 2.15. 

- S,„) = - S,„) (2.17) 

The above values result in Obukhov lengths of 0.59 m for the maximum friction 
velocity or 0.20 m for the average friction velocity. Assuming the fresh water 
uniformly enters the ocean (by using the spatial average of surface melt rate of 1.5 
cm day'*) the salt flux is reduced to a lower bound of 4.6x10'^ kg m'^ s'*. In this 
case, Obukhov length was 8.9 m or 3.1 m for maximum or mean surface stress. 

The Obukhov length estimates indicate that when the meltwater drained 
laterally into leads (upper bound) in the early summer, turbulent stress was too 
weak to mix it vertically: Obukhov lengths were less than the ice draft. Even 
using the average surface stress over days 190-207 and horizontally averaged 
fresh water flux (the lower bound) results in an Obukhov length (3 m) on the 
order of ridge keel depths. Only during periods of above-average surface stress 
and away from leads did the Obukhov length exceed the nearby ridge keel drafts. 
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The ratio of the neutral-stratification eddy size {Q.Q5u*o/f) to the Obu kh ov 
length is the “turbulent lead number” of Morison et al. [1992], The turbulent lead 
number was developed for use in unstable convection from wintertime leads, but 
its basis in the turbulent kinetic energy equation does not require a particular sign 
for buoyancy flux. For freshwater fluxes in the lead near the upper bound (7x10'^ 
kg m' s' ), this lead number is less than one, implying that buoyant destruction of 
turbulence dominated the entire boundary layer in the lead. In other words, 
turbulent eddies generated by ice motion were inhibited by the stabilizing effect 
of meltwater at the lead surface. AMTV observations from day 204 (Fig. 2.8) 
support the low turbulence levels below 4 m in the lead: at 6 m the vertical water 
velocity and temperature fluctuations were small and heat flux at 6 m was very 
small or zero, despite the solar heating of the lead surface. Friction velocities 
were among the smallest observed with the AMTV in the mixed layer. 

The persistence of the fresh layer in the lead can also be examined in 
terms of the internal Froude number: 

Fr = . ^ where 77/ is the layer thickness and g = ^ (2.18) 

^g'H, P 

We assume the layer moves with the ice and is contained by the ice draft. (The 
draft exceeds the layer depth.) When the velocity of the fluid below the layer 
(approximately the ice velocity) exceeds the wave speed at the interface, the ice 
can no longer block the flow of fresh water downstream. Under this supercritical 
condition, the layer begins to slide under the ice, is forced past a rough surface, 
and is entrained into the mixed layer. Prior to day 207, ice speeds were less than 
10 cm s'' and the layer salinity was less than 5 psu, compared to 31 psu for the 
seawater. For the Froude number to be greater than one and allow the layer to 
escape, the layer would have to be less than 0.06 m. The lead layer probably 
reached this depth or greater during a few hours of low ice velocity and survived 
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the next significant ice movement. In terms of Fronde number criticality, as the 
layer grew it became more resilient to ice motion. By day 185, the interface was 
at a depth of 20-30 cm \Paulson and Pegau, 2001]. On day 208, the ice velocity 
finally exceeded 10 cm s'*. By this time, the layer was 1.2 m thick. In this case 
Fr= 0.3 (ice speed of 0.12 m s'*, layer thickness 1 m, and salinity jump of 26 psu). 
The Froude number suggests the fresh water in the lead was still stable, and 
indeed the layer remained in the lead as the ice began to move (until day 212.5 
according to the ROV data). 

However, the layer was shrinking at a rate of about 13 cm day'* from day 
206.5 to day 210.5 (based on ROV observations). The lead-edge TIC observed 
large heat and salt fluxes consistent with warm, fresh water mixing downwards on 
day 208 (Figs. 2.17 and 2.19). Apparently other mechanisms were acting to thin 
the fresh layer while the ice velocity remained sub-critical. It is likely that as the 
ice began to move on day 206.5, the layer was eroded away by turbulence 
generated upstream of the lead more quickly than the melt was supplied. After 
day 210 ice divergence also thinned the layer: ice divergence in a 50 km by 50 
km box containing SHEBA became positive on day 210 (2.5% per day) and 
remained positive for the next several days [Stem and Moritz, 2002]. On day 
211, the thickness was 36 cm, and the ice speed was up to at least 0.20 m s'*, so 
that F?^0.8. By the end of the day the layer was at 12 cm, and it was gone the 
next day. The order-one Froude number on day 211 suggests the layer finally 
became thin enough to be unstable and was rapidly flushed out of the lead. From 
day 206.5 to day 212.5, the fresh lead layer was destroyed by a combination of 
erosion from below, ice divergence, and hydrodynamic instability. 

2.4.2 PERIOD II: MELTWATER ENTERS OCEAN 

The transition stage at SHEBA (210-214) finally saw the mixed layer 
freshen dramatically (Fig 2.5a). The observed mixed layer change implies a 
surface salt flux over that period of 4.6x10'^ kg m'^ s'*. To calculate the salt flux 
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from mixed layer salinity, one must assume that the change in salinity at 30 m is 
unaffected by the surface and that in the absence of surface processes the salinity 
profde above 30 m would hold its initial shape throughout the period. The latter 
assumption is questionable in such a variable region as the Chukchi Cap, i.e. there 
could be advective effects above 30 m that are not identical to those at 30 m. On 
the other hand, the total amount of freshening observed at SHEBA (minus the 30 
m advective change) roughly corresponds to the total ice and snowmelt: 112 cm 
of fresh water similar to the observed 126 cm. Perovich et al. [2003] report the 
latter value for days 147-244. In addition, buoys about 20 km from SHEBA 
showed simultaneous mixed layer freshening of the same magnitude as the 
advection-corrected SHEBA CTD [McPhee, personal communication]. The 
freshening occurred 2-3 days before the SHEBA camp left shallow water, which 
also suggests that much of the mixed layer freshening was not linked to the drift 
into deeper water. 

The flux implied by mixed layer freshening is much larger than the heat 
balance estimate (Fig. 2.4) and the lead flushing flux. From section 2.3.2, the lead 
salt fluxes were 5x10'^ kg m'^s'' for day 206.5 to day 210.5 and 12x10'^ kg m'^ s’' 
for day 210.5 to day 212.5. Since the lead fraction was about 5%, the mixed layer 
would experience fluxes about 5/95 or 5% of these magnitudes, not nearly enough 
to account for the observed salinity changes. The heat balance estimate of salt 
flux is also too small and shows little correspondence to changes in the mixed 
layer. Clearly surface retention and subsequent runoff were a major factor 
[Eicken et al, 2002]. Using 4.6x10' kg m' s' and the maximum friction speed 
of 0.015 m s'^ results in an Obukhov length for the storm period of 24 m. Using 
an average friction speed of 0.01 m s'^ results in a value of 7.1 m. As expected, 
the high levels of shear-generated turbulence during the storm reached the mixed 
layer and deepened it despite the very large stabilizing buoyancy flux (Fig. 2.5b). 

During the high winds on day 213 the AMTV showed particularly strong 
friction speed at 8 m (Fig. 2.11): the run average of 0.012 m s'* was larger than 
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the surface value (Fig. 2.16). Heat flux was also much higher than the bulk 
estimate, averaging 60 W m' upward over the run (Fig. 2.17). The most 
interesting result is the strong downward salt flux, averaging -1x10'^ kg m"^ s'^ 
(Fig. 2.19), despite a flux of +4.6x10'^ kg m'^ s'^ estimated from ocean freshening 
during days 210-215. Most of the activity took place under and downstream of a 
ridge (Fig. 2.11). A second run also at 8 m a few hours later showed even 
stronger upward heat flux and downward salt flux under regions of ridges and 
small leads (run averages of 0.015 m s'\ 100 W m'^ and -2x10'^ kg m'^ s'*). 
Because these two runs took place when the lead fresh layer had just disappeared 
(212.6 and 212.8), it is possible that the fresh warm plume from the lead was 
forced under adjacent ridges and under saltier, cooler water from the downstream 
side of the ridge. In general, this is a mechanism that may arise when a strong 
downward fresh water flux is upstream of a weak or small fresh water flux. If the 
less dense fresh water from the lead is forced downward as it is swept under the 
adjacent ice, it is possible that it undercuts denser, saltier water. This unstable 
situation could only occur if the horizontal density gradient has the same sign as 
the horizontal fluid flow relative to the ice, e.g. at the downstream edge of a 
summer lead. Such “over-running overturns” would enhance the turbulent mixing 
at the downstream lead edges causing an increase in the variance of vertical water 
velocity (implying a higher ice-ocean stress and possibly a deeper mixed layer) 
and an increase in its characteristic length scale. See Crawford et al. [1999] for a 
discussion of observations suggesting the presence of this mechanism in the 
mixed layer under land-fast ice. The over-running overturns may explain the 
downward salt flux and strong heat and momentum fluxes observed by the 
AMTV in this region. 

2.4.3 PERIOD III: QUASI-STEADY BOUNDARY LAYER 

The period after day 215 was characterized by weak freshwater fluxes and 
stronger winds. A brief period of low wind stress (around day 215) may have 
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allowed the surface melt to build up slightly in a layer in the lead, but not to a 
thickness to which it could be maintained through the next wind event. On day 
215, the friction velocity and salt flux at the surface were 0.005 m s'^ (Fig. 2.16) 
and 4x10'^ kg m'^ s'* (Fig. 2.19), respectively. The surface salt flux should be 
roughly a factor of two larger than the AMTV flux at 5 m in this range, according 
to similarity theory (Eqs. 2.13-2.15). The resulting Obukhov length for day 215 is 
about 1 m. By day 220, the friction velocity had increased to nearly 0.01 m s'*, 
and the salt flux at the surface was about 2x10'^ kg m'^ s'*. The Obukhov length 
increased to 16 m. No fresh layer was observed. 

Despite highly permeable ice, the meltwater found the lead surface before 
mixing downward, presumably due to surface and internal drainage. Indeed, on 
day 220 the AMTV still observed the salt flux to be in the lead, not under the ice 
(Fig. 2.14). The heat flux at 5 m in the lead was downward. The lead surface was 
being warmed and freshened, and these properties were mixed downward by 
turbulence advected from upstream and generated by wind at the lead surface. 
The cumulative effect of the vertical convergence of fluxes on a particular water 
parcel in the lead was to warm and freshen it as it neared the lead edge (Fig. 2.12). 
At downstream distances greater than about 100 m under the ice, fluxes were 
smaller or did not pass the significance test, and the friction speed at 5 m 
increased. The variances of temperature and salinity were also reduced in the 
region (end of segment 3 and beginning of 4). The data seem to describe a steady 
internal boundary layer whose structure was determined from the time dependent 
surface boundary conditions of heat and salt fluxes, and stress. 

One can make a simple estimate of the depth, h, of the internal boundary 
layer (IBL) as a function of distance downstream, x, from the change in roughness 
length. It is based on the change in roughness length from the lead “surface” to 
the ice surface [Garratt, 1990; Kantha and Clayson, 2000]. 
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h , ,h ^ h , r. 1 X 

—ln(—) + \ = 2k^ — 


(2.19) 


Simply put, one assumes that the “signal” of the new boundary roughness 
“diffuses” via turbulence away from the boundary with distance downstream, 
much like a spot of dye injected at the wall of a boundary layer flow would 
diffuse. In our case the signal is momentum. The bulk of the IBL work has been 
applied to the neutral atmospheric boundary layer, in which case one can safely 
assume as in the references cited above that the IBL is within a surface layer. In 
the ice-ocean boundary layer, the surface layer is much thinner, certainly less than 
5 m, so the solution above is not appropriate once the IBL reaches this thickness. 
If one assumes the IBL continued to grow at the maximum rate found at the outer 
edge of the surface layer, one can constrain the IBL depth for a given distance 
downstream. For day 220, assuming a roughness length of 0.01 m, the IBL 
reaches a depth of 5 m at a distance of at least 81m downstream, close to the 
visual estimate of 100 m in Fig. 2.12. The agreement with this simple scaling 
strengthens the argument for vertical mixing of momentum. It seems that the 
boundary layer at this time can be approximated by a one-dimensional 
(horizontally homogeneous), unsteady boundary layer. 

2.4.4 BULK, TIC, AND AMTV INTERCOMPARISON 

As described, the bulk estimates of surface fluxes disagree substantially 
with the TIC data, even when adjusted for attenuation with depth. See section 
2.3.3 and Fig. 2.15. The attenuation factor was derived through a similarity 
solution for flux profiles in a boundary layer stabilized by buoyancy flux. The 
heat balance estimate of salt flux was used, which was seen to be too small to 
account for mixed layer salinity changes. Assuming a surface buoyancy flux that 
includes an estimate of runoff more strongly attenuates the surface friction speed 
and heat flux and produces good agreement between day 211 and day 215, 
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especially at the lower cluster (not shown). However, this buoyancy flux is an 
order of magnitude larger than any data in Fig. 2.15c. The TIC salt flux is called 
into serious question since it is far less than the flux required for the observed 
salinity changes and is inconsistent with the TIC-observed heat flux and stress. It 
is possible that biofouling led to poor conductivity measurements and/or that the 
turbulent scales were too small to be resolved by the SBE-04 conductivity cell. 

On day 219 and day 220, it appears the TIC data were more reliable. 
Average surface buoyancy flux was probably much smaller as the mixed layer 
salinity was fairly constant over the period following day 215. It could be that the 
TIC being cleaned more regularly or the lack of biological matter combined with 
the smaller buoyancy flux (and so larger scales) allowed the conductivity cell to 
make an accurate measurement. If so, the fresh water flux from surface melt was 
enough to cause the actual flux to significantly exceed the heat balance estimate. 

The AMTV-derived heat flux qualitatively makes sense. Before day 215, 
most of the AMTV data segments were obtained under ice. Very few runs 
showed significant heat flux before the wind increase on day 207. Between day 
207 and 215, the AMTV generally showed mixed layer heat flux values less than 
the surface value but in the same direction (upward) consistent with the 
attenuation of the heat flux with depth predicted by similarity theory. The 
exception is day 213, see the discussion of over-running overturns in section 
2.4.2. After day 215, when most of the AMTV data were obtained under open 
water, heat fluxes in the mixed layer were downward because of radiative flux at 
the surface and the lack of ice to absorb heat by melting. Of course the ocean city 
TIC did not indicate downward fluxes as the site was far from open water. 

Finally, we note that the AMTV-derived friction velocity agrees very well 
with the steady-state force balance estimate. The AMTV samples many surface 
types, so it provides a more representative measurement of spatial-average stress 
than a point measurement. Point measurements of stress close to the underside 
of sea ice can be difficult to use as representative of large-scale values 
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[Wettlaufer, 1991; Morison et al., 1987], As the lead-edge TIC data confirm, the 
choice of location for a point measurement of stress may result in different values 
for friction velocity. At SHEBA the lead-edge stresses were not representative of 
the mean conditions. It is important that in a spatially varying boundary layer the 
AMTV can accurately measure ice-ocean stress in addition to exchanges of heat 
and salt. 

2.4.5 NUMERICAL MODEL MOTIVATION 

The first-order closure model of McPhee [1987; 1992; McPhee and 
Kantha, 1989] is described below and used to make the case that vertical turbulent 
“diffusion” driven by a time-varying surface condition dominate the quasi-steady 
period. The model will be used to judge whether the observed heat and salt fluxes 
and characteristics of vertical water velocity (variance and characteristic 
wavenumber) are quantitatively consistent with the one-dimensional view. The 
model is also used to investigate the largest oceanic flux event of the su mm er 
from day 207-215. As we will see, observations at SHEBA on days 207-215 
indicate the observed vertical mixing was deeper than predicted by the one¬ 
dimensional model. During the period, the large salt and heat fluxes were likely 
to be the strongest under leads. If there were a time when horizontal gradients 
were important, this would be it. The three-dimensional time-dependent model 
described below is required to investigate the importance of the new mixing 
mechanism in both the strong and weak fresh water flux scenarios. 

We begin with an analytical description of the ice-ocean system that 
expresses conservation of momentum, temperature, and salinity. A description of 
a 1-D model, its previous uses, and results of boundary layer simulations follows. 
A 3-D model is described and compared with boundary layer observations. 




remotely-operated vehicle (ROV—tethered), conductivity-temperature-depth (CTD) instruments, and turbulence instrument 
clusters (TIC). Other observations were: solar insolation, wind and ice velocity, aerial photography, and CTD and TIC 
deployments away from the lead. 
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Figure 2.2. Drift of the SHEBA ice camp from October 2, 1997 to September 26, 1998. June, July, and August are blue, cyan, 
and black segments over Chukchi Cap. Positions: Des Groseilliers logbooks. Bathymetry: http://www7430.nrlssc.navv.mil . 
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Figure 2.3. Outline of the SHEBA lead known as “Sarah’s Lake” on three days in 
1998. The location of the operations hut is at the origin. Lead perimeter was 
collected by tracking an acoustic beacon with the Applied Physics Lab tracking 
range. The beacon was towed with a small boat. The lead on day 217 was not 
mapped in its entirety. 
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Figure 2.4. (a) Friction speed (blue) at the ice-ocean interface from the SHEBA 
Project Office wind speed at 10 m (red) and ice velocity, assuming a steady state 
balance, (b) Elevation above freezing temperature at 4 m depth, (c) Ice-ocean 
heat flux from bulk exchange relation, proportional to the product of (a) and (b). 
(d) Salt flux at ice-ocean interface assuming the bulk heat flux goes to melting ice 
of salinity 2 psu. 
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Figure 2.5. Ocean city CTD profiler salinity and ocean depth (a), buoyancy 
frequency in top 40 m (b), potential temperature (c), and potential temperature 
departure from surface freezing temperature (d). In (b), the 1x10""^ s"^ contour is 
the thick magenta line, and the black circles are AMTV run depths and times. 
































46 



0 




(c) 

/ 

i 


- 1.7 - 1.6 - 1.5 - 1.4 

Temperature (“C) 



□ 






■•(b).... 



5 10 15 20 25 30 35 

Salinity (psu) 



30 30.05 30.1 30.15 30.2 

Salinity (psu) 


Figure 2.6. Remotely-Operated Vehicle profiles of temperature (a) and salinity 
(b) in the main SHEBA lead. On 7/25/98 (day 206) and 7/29/98 (day 210) the 
ROV could not float passively to the surface due to the strong stratification in the 
lead, so the profiles do not extend to the surface. The profiles on 8/7/98 (day 219) 
are shown in detail in (c) and (d), along with surface freezing temperature. 
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Figure 2.7. AMTV run on day 204.0 (UT). (a) Vehicle depth marked with 
symbols representing the start and finish of each leg. Also shown is the ice draft 
observed by the AMT Vs upward-looking acoustic altimeter, (b) Temperature 
observed by the AMTV. (c) Horizontal position of the AMTV as recorded by an 
acoustic tracking range. The x andy axes are defined in the tracking range 
software and remained constant throughout the experiment. The same symbols as 
in (a) are used to indicate the segments of the run to be analyzed. The outline of 
the lead is shown, as well as the ice velocity of 6.4 cm s '. The open water returns 
from the altimeter are plotted with thicker line segments. 
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Figure 2.8. AMTV run on day 204.0, turbulent quantities, (a) Detrended 
temperature for each leg and ice draft, (b) Vertical water velocity calculated from 
vehicle motion, (c) Heat flux with 100 m bin averages shown as horizontal bars. 
Only the two indicated are significant (10 and -3 W m'^). (d) Friction speeds from 
the spectra of vertical water velocity and the 95% confidence interval. 
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Figure 2.9. AMTV run on day 213.0 (UT), as in Fig. 2.7. At this time the lead 
was nearly closed. As before, the solid arrow indicates the ice velocity. The 
AMTV crossed a ridge at around +200 m on the x-axis. 
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Figure 2.11. AMTV run on day 213.0, as in Fig. 2.8, except salt flux is also 
plotted. Significant heat flux is observed at the 500 m, 600 m, 700 m, and 900 m 
bins (73, 72, 31, 74 W m'^). Significant salt flux is observed at the 500 m, 600 m, 
900 m, and 1000 m bins (-1.1, -1.3, -1-0.82, -1.0x10'^ kg m'^ s'*). These segments 
are under and downstream of the ridge as indicated by the dashed line. 
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Figure 2.14. Fluxes for AMTV run on day 219.9. As in Fig. 2.11. Most of the 
bins show statistically significant fluxes as indicated. The lead fluxes of heat and 
salt averaged 98 and 144 W m'^ and 1.7 and 2.4x10'^ kg m'^ s'* (segments one and 
two). Segment three did not pass the significance test for heat or salt flux, and 
segment four only had significant salt flux: 1.9x10'^ kg m'^ s'*. The segment 
averages reported here are calculated from the covariances over each entire 
segment, not the average of the 100 m segments. The first 200 m of segment one 
were not included as it appears to have been contaminated by salinity drift (see 
Fig. 2.13). 
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Figure 2.15. (a) Turbulence Instrument Cluster (TIC) observed friction speed at 
3.5 m (red) and 7.5 m (green) compared to surface friction speed from force 
balance relation. The small dots are friction speed at the cluster depths 
extrapolated from the surface value using similarity theory, (b) TIC observed 
heat flux, compared to bulk estimate of surface heat flux and similarity 
predictions, (c) TIC observed salt flux comparison. Surface value based purely 
on thermal balance (no runoff). 
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Day of 1998 


Figure 2.16. Ice-ocean interface friction speed from steady-state balance of wind 
stress and Coriolis force compared to all AMTV run averages (pink crosses), 
AMTV run averages using only boundary layer points (red squares), and lead- 
edge TIC observations (3.5 m blue stars and 5.5 m green stars). 
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Figure 2.17. Ice-ocean heat flux from bulk estimate, and heat flux in the mixed 
layer from AMTV, and lead-edge TIC. Bulk heat flux is based on elevation of 
mixed layer above freezing and surface friction velocity from Fig. 2.16. Only 
values of AMTV heat flux that passed the significance test are shown. After day 
215, the AMTV was predominantly in open water where downward heat fluxes 
were observed. 
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Day of 1998 

Figure 2.18. Detailed view of heat flux on days 219-220. Blue line with dots is 
the downward shortwave radiation measured just above the ice surface by the 
SHEBA Project Office. AMTV and lead-edge TIC heat fluxes correspond well to 
the radiative flux. 
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Figure 2.19. Salt flux in the upper ocean at SHEBA. Bulk estimate (blue line) is 
the melting due to the bulk heat flux in Fig. 2.17. The conductivity cell was not 
operating correctly before day 213. Only AMTV values that passed the 
significance test are included in the averages shown. The AMTV sampled mostly 
open water after day 215 where upward salt flux was typically observed by both 
instruments. 
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3. MODELING STUDIES 

3.1 ANALYTICAL FORMULATION 

Both models in this study use conservation equations written for the mean 
quantities with the vertical coordinate positive upward. An incompressible, 
Boussinesq fluid is assumed. Reynold’s decomposition is carried out, which 
means that there must be a separation in scale between a “mean” signal and a 
“turbulenf’ signal (the latter indicated by a single prime). The momentum, 
temperature, and salt fluxes are expressed as covariances among the turbulent 
variables. The method used here to determine these covariances from mean 
properties (turbulent closure) is first order closure, meaning the fluxes are related 
to the gradients in the mean properties through turbulent diffusion coefficients. 
The primes on pressure and density indicate deviations from the mean, stationary 
state. The mean vertical density profile and associated pressure profile are 
subtracted from their actual profiles. 


Ou - Vp' 

— + u-Vu + ykxu = —— + V • T - 
dt 


A 


(3.1) 


where t, =-{u'.u'j)nj 


is the /-component of stress and n is the unit vector 


(summation over j is understood) 


V-u = 0 


(3.2) 


— + u-VT = -V-¥^ + Q^ (3.3) 

dt 

where Ft is the kinematic heat flux: Fy, = (wT')i + (vT')j + (wT')k 

/ z 

and Qt is a radiative heating term: = —exp(—), where L is an extinction 


depth equal to 4 m and 4 is the shortwave flux at the ocean surface. 
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— + u-VS = -V-¥, (3.4) 

at 

where Fs is the salt flux: = (u'S')i + (v'-S")] + 


Boundary conditions depend upon the particular problem of interest, but in 
general, we will specify fluxes of momentum {u*o), heat and salt 

at the ice-ocean interface or the atmosphere-ocean interface (for open 

water). 

In general, density depends on temperature, salinity, and pressure through 
the nonlinear equation of state, but for our purposes density is a simple linear 
function of salinity: 

p = 0.8085+ 1000 (3.5) 

It is acceptable here to neglect effects of temperature on density because near the 
freezing point, the thermal expansion coefficient of seawater is much smaller than 
the saline contraction coefficient. Pressure effects are small because we are 
concerned with surface mixed layer processes. 

3.2 ONE-DIMENSIONAL MODEL 
3.2.1 MODEL DESCRIPTION 

Results from previous studies show that in many cases, a one-dimensional, 
first-order closure model for the stratified oceanic boundary layer simulates 
observations well. In particular, McPhee [1987] and Mcphee et al. [1987] 
successfully used similarity scaling arguments and a time-dependent 1-D model to 
describe the under-ice boundary layer during the 1984 Marginal Ice Zone 
Experiment (MIZEX). His model [McPhee, 1999] agreed well with a subset of 
the observations and with the results of the level 2.5 model of Mellor and Yamada 
[1974; 1982]. McPhee and Kantha [1989] extended the model to include internal 
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wave generation by sea ice, which helped explain the high ice-ocean drag during a 
later period of MIZEX. Morison et al. [1992] used the model to examine unstable 
“forced” convection from leads, with reasonable results. McPhee [1994, 1999] 
demonstrated some interesting features of Weddell Sea mixed layer development 
with the one-dimensional model. The model of this section is similar to that of 
McPhee [1987] and is recapitulated here with some modifications. 

The model assumes purely horizontal and horizontally homogeneous flow. 
Velocity is expressed in a reference frame moving at a far-field velocity (beyond 
influence of the boundary). The vertical velocity vanishes because of Eq. 3.2 and 
a w=0 top boundary condition. The mean momentum equation then reduces to 
the following. 

— + ykxu = — = — (K —) where u is the horizontal velocity (m,v,0) (3.6) 

dt dz dz dz 

The conservation equations for temperature and salinity are: 


dt 


d 

dz 


{w'r) + Q, 


d , „ 

{oC'fK 

dz 


dz 


) + Qt 


dt 



d , „ 

dz 



(3.7) 

(3.8) 


Vertical fluxes are modeled as vertical property gradients multiplied by an 
eddy exchange coefficient, K. For temperature and salinity, this is multiplied by 
the ratio of turbulent heat or salt diffusivity to eddy viscosity (the turbulent 
Prandtl and Schmidt numbers, respectively). A first-order turbulent closure 
model [McPhee, 1981; 1994] specifies the eddy viscosity as the product of 
characteristic turbulent velocity and local turbulent mixing length: 


K = u„A, 


(3.9) 



63 


The local friction velocity, u*, is the square root of the local stress, calculated 
from a previous estimate of eddy viscosity and the velocity shear profile. 


u» = 




(3.10) 


The friction velocity is set to a constant value in the surface layer. 

McPhee [1994] develops the mixing length calculation for the stable or 
neutral mixed layer. In the inner layer, observations support the assumption of 
linear growth of the mixing length, X, with distance to a maximum value 
[McPhee, 1981]. The maximum value in the mixed layer is calculated based on 
interfacial surface stress and planetary rotation and is bounded by stabilizing 
buoyancy flux through the stability parameter, ;/*, which is less than or equal to 
one. 

T = min(/f|z|, (3-11) 


rjl — (1 + 


xfR-cLo 


(3.12) 


The mixing length retains the constant value throughout the mixed layer, below 
which it is reduced in the strong stratification of the pycnocline. There the mixing 
length is calculated from Eqs. 3.11 and 3.12 using the local buoyancy flux and 
local stress. The closure expressions are arrived at through scaling the steady 
momentum equation appropriately for the boundary layer, making use of an 
assumption that the nondimensional eddy viscosity is a universal constant away 
from the interface (in the outer layer) [McPhee, 1981]. 
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Turbulent Schmidt and Prandtl numbers are calculated using a gradient 
Richardson number approach, outlined in McPhee [1994], They are assumed to 
be equal in this model, and when vertical shear is strong, they are equal to one: 


OC g — OCf 


X R/< 0.079 

^ exp(-1.5VRi-0.079), 0.079 < Pi < 5 
0.039, Pi > 5 


(3.13) 


g 

A 

dz 


5u 

2 


dz 



(3.14) 


The boundary layer is viewed as a single column, with surface boundary 
conditions varying as the ice-lead system passes over. Boundary conditions at the 
surface are heat flux, salt flux, momentum flux (stress), and roughness, which can 
all vary in time. Stress is specified in the model in the form of surface friction 
speed, u*o, and roughness length, zo, is also specified. Heat flux at the interface is 
made up of two terms: the downwelling shortwave radiation and ice-ocean heat 
flux when water temperatures above freezing in the mixed layer contact the ice 
boundary. The former is distributed in an absorbing layer via the heating term 
with an extinction depth of 4 m. The latter is calculated from the bulk heat flux 
transfer model (Eq. 2.8) [McPkee, 1992]. Salt flux is also made up of two terms: 
a percolation velocity and salt flux due to interfacial melting. The former is a 
specified velocity of the interface (upward) due to fresh meltwater entering the 
top via runoff or percolation through the pack ice. The latter component of salt 
flux is calculated using the salinity boundary relation of Eq. 3.15 (a more 
complete version of Eq. 2.9) McPhee [1990]. Conduction through the ice is 
assumed negligible here. 
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{^^fS\^{w'r\ASq;' +F„.,,,-') (3.15) 

where AS = S^,-Sj^^ and =83.8(1-0.035',.^g). Bottom boundary conditions 

are zero flux conditions for momentum, temperature, and salinity. 

Initial conditions on temperature and salinity profiles are set, usually 
according to an observed set of profiles. Initial conditions for velocity and eddy 
viscosity are not as straightforward. The steady version of the one-dimensional 
momentum equation is solved with an initial estimate of eddy viscosity. This 
solution is then used to recalculate mixing length and friction velocity profiles, 
together giving a new eddy viscosity profile. The steady solution is calculated 
again, and this process is iterated to convergence for initialization of velocity and 
eddy viscosity. 

It is possible to view the 1-D, horizontally homogeneous unsteady 
boundary layer as a 2-D steady boundary layer. Mellor et al. [1986] and Morison 
et al. [1992] apply the transformation discussed here. The boundary conditions 
and boundary layer equations are transformed to a coordinate system moving with 
the ice in the negative x direction at constant speed t/,. Let x* represent the 
horizontal distance in the x direction downstream of an arbitrary point in the ice 
pack: 

x* = C/.t + x (3.16) 

Equation 3.6 becomes 

^,^ + /kxu = —(i:— ) where u = (m,v,0) (3.17) 

ox * oz oz 


It is important to realize that Eq. 3.17 does not violate the horizontal homogeneity 
assumption because it is only a transformation of Eq. 3.6. A vertical profile at 
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some value of x* is still assumed to be the solution to the horizontally 
homogeneous equations. It depends on the boundary conditions at that location 
and the solution at a neighboring value of x* that acts like the solution at the 
previous time step. It represents the spatial variation associated only with the 
temporal variation of the 1-D boundary layer profde. Gradients in the new 
horizontal variable are not used in the advective terms of the momentum equation, 
so for an accurate solution these gradients must be small. Unlike time-dependent 
solutions, those expressed in terms of x* can be compared directly to AMTV data, 
since they are both expressed in horizontal distance related to the ice pack. 

3.2.2 RESULTS AND DISCUSSION—DAY 220 

A two-dimensional steady simulation of the lead around solar noon of day 
219, 1998 is shown in Fig. 3.1. The model domain is 4.25 km in the horizontal by 
30 m in the vertical with 0.15 m spacing in the vertical and 5 min time steps 
(corresponding to 51 m in the horizontal for the average ice velocity of 0.17 m s'* 
to the left). The model lead is a 1.25 km segment in the upper left with incoming 
shortwave radiation of 448 W m' , extinction depth of 4 m, and albedo of 0.08. 
The lead is further specified by a surface friction velocity of 8x10'^ m s'* towards 
the left, a roughness length of 2.7x10'^ m, and a fresh water input equivalent to 
0.5 cm day * surface melt over an 80% ice-covered ocean draining into the 
remaining 20% lead fraction (a salt flux in the lead of bxlO'*’ kg m'^ s'*). The 
surface stress and roughness length were calculated using the observed 10-m wind 
speed and techniques from Steele et al. [1989]. The fresh water input due to 
surface melt is determined from measurements of the ablation rate [Perovich et 
al, 2003]. The salinity of the ice is taken to be 0.5 psu [Eicken et al, 2002]. The 
top boundary conditions under ice downstream of the lead do not include any 
source of meltwater other than bottom melt, and the radiative flux is zero. The 
surface friction velocity and roughness length increase to 0.01 m s'* (still towards 
the left) and 0.01 m, calculated from observed ice velocity using a Rossby 
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similarity law and a representative value for sea iee, respeetively. The initial 
temperature and salinity profdes represent those measured at the lead edge during 
the AMTV run with the Sea-Bird SBE-19 CTD: a 15-meter mixed layer with 
temperature of more than 0.1 C° above the freezing temperature. 

Results show that, like the AMTV observations, the temperature at 5 m 
inereases steadily as the downstream lead edge is approaehed from upstream, then 
falls and levels off under the ice. As discussed above, this is consistent with 
continuous solar heating and mixing down of the surface water. Figure 2.12 
shows the temperature observed (also around solar noon) by the AMTV at 5 m 
increased at a rate of about 0.01 deg per 100 m in the lead, while the model results 
indicate a rate of about half that value. Also, the 5-m observed temperature 
decreases downstream of the ice edge more quickly: 0.02C° over 100 m 
compared 0.02C° over 1000 m in model. Despite the difference in horizontal 
temperature gradient, the model reproduces a qualitatively correct temperature 
structure. 

The model heat flux in the lead at 5 m is about 100 W m'^ downward over 
the 500 m just before the lead edge, compared to an average of 98 W m' and 144 
W m' for the AMTV over the two large lead segments on day 219.9. 
Downstream, under the ice, the model heat flux at 5 m strengthens as the rough 
ice surface rapidly mixes down the warmer surface water. The heat flux at 5 m 
then relaxes slowly as the vertical temperature gradient weakens, eventually 
becoming upwards about 2 km downstream of the lead edge. The heat flux at 10 
m remains downward throughout the run and is more downward than the 5 m heat 
flux between 1 km and 2 km. The AMTV heat flux observations indicate a much 
faster decay of heat flux under the ice. While the 100 m of AMTV data just 
downstream of the edge show a strong downward flux of 97 W m' at the very end 
of the run, the segment 100-400 m downstream of the lead edge averaged 40 W 
m'^ and there were two segments with insignificant fluxes. The model downward 
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flux remains greater than this level until about 1.2 km downstream of the lead 
edge. 

The salinity at 5 m shows a freshening also consistent with vertical mixing 
of fresh water input at the surface of the lead. The freshening trend weakens 
under the ice, since the source of fresh water (basal melting) is weaker than the 
surface melt drainage flux in the lead. The absolute salinity observations from the 
AMTV in Fig. 2.12 have been contaminated by a low frequency drift, which 
makes comparison with the model difficult. However, the trend in salinity from 
the turn-around point at 800 m is in the same sense as the model: freshening as 
the lead edge is approached from upstream. Any trend before the turn-around 
point was lost in the attempt to correct for the sensor drift. The salinity deviations 
are unaffected by the low-frequency sensor drift, allowing the calculation of salt 
flux. The salt flux in the last 500 m of the model lead averages to 0.4x10'^ kg s'' 
m'^ upward. The AMTV salt flux was 1.7x10'^ kg s''m'^ and 2.4x10'^ kg s''m'^ 
upward for the two 500 m lead segments. The model seems biased low; perhaps 
the freshwater flux in the lead was greater than the value used in the lead 
boundary condition. In the 100 m downstream of the lead edge, the model salt 
flux under the ice at 5 m increases to about 0.8x10' kg s' m' and decreases 
gradually. The observed under-ice flux on day 219.9 (taking the covariance over 
all of segment 4) gives 1.9x10' kg s' m' , while the segment 3 covariance does 
not pass significance test. Besides the difference in flux magnitude between the 
model and observations, the model shows an increase in salt flux relative to the 
lead, while the observations indicate a moderate decrease relative to the lead. 

The modeled mixing length and stress can also be compared to the 
observed. The modeled mixing length at 5 m is about 0.5 m in the lead, 
increasing to 0.8 m under the ice and decreasing gradually with downstream 
distance. The increase in modeled length scale from lead to ice boundary 
conditions is due to the reduction in stabilizing buoyancy flux and increase in 
stress at the surface. The observations do not show a clear pattern, but they are in 
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rough agreement with the model. By Eq. 2.4, the mixing length in the lead at 5 m 
is 0.7 m with a maximum under the ice edge of 1.8 m. 

The model shows an increase in friction speed at 5 m from 0.6 cm s'^ to 
0.8 cm s'\ corresponding to the specified jump in surface friction from open 
water to ice of 0.8 to 1.0 cm s'*. This is in good agreement with AMTV results. 
The AMTV friction speeds on day 219.9 calculated from the vertical water 
velocity spectra and shown in Fig. 2.14 were around 0.8 cm s'* in the lead and 1.0 
cm s'* under the ice. 

Contours of model density are shown in Fig. 3.2. The lead is in the upper 
left and the surface stress is to the left. The contour of zero horizontal density 
gradient is overlaid, roughly separating the domain into one region influenced by 
the lead surface and another by the ice surface. The horizontal density gradient is 
positive in the region under the ice, a situation conducive to the over-running 
overturns discussed previously. However, the 2-D steady model does not allow 
these overturns, since the advective terms have been neglected. Even in the 
quasi-steady period around day 219, it appears that there exists a potential for 
instability that cannot be simulated with the current model. 

The model represents the boundary observed on day 220 reasonably well. 
Both the model and AMTV observations indicate fresh, warm water mixing 
vertically in the lead. However, the observations suggest the temperature, 
salinity, and heat and salt fluxes diminish significantly within a few hundred 
meters of the downstream lead edge, while the model values decay over much 
larger distances. One possible reason for the differences in horizontal dependence 
of both scalars and fluxes is that the advective terms are in reality significant, 
particularly near the lead edge, causing over-running overturns and enhanced 
mixing. However, unstable conditions were not directly indicated by the day 220 
AMTV observations. That may be because the horizontal density gradients were 
too weak to play a significant role. The enhanced vertical mixing observed under 
the ice may also be because the model boundary condition (stress and roughness 
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length at 50 m intervals) could not adequately capture the effect of the lead edge 
on mechanical mixing. Lead edge mechanical forcing was suggested to be an 
important factor in the pattern of unstable lead convection by Smith and Morison 
[1998], The 2-D unsteady model is required to investigate these two possibilities. 

3.2.3 RESULTS AND DISCUSSION—DAY 209 TO DAY 215 

While the one-dimensional model produced reasonable simulations for the 
boundary layer transition from lead to ice on day 220, we now investigate its 
ability to simulate the horizontally averaged boundary layer during the strong 
fresh water flushing period from day 207 to day 215. There is no lead in this 
simulation. The model domain is 30 m deep with a vertical grid spacing of 0.3 m. 
The time step is 15 minutes. Initial temperature and salinity profiles for the 
model are from the cast at ocean city on day 209.2. Initial velocity and eddy 
viscosity profiles are calculated using the steady solution for the initial boundary 
condition. The model is forced by surface stress from Fig.2.4a. Radiative flux 
through the ice is set to zero, and the ice-ocean heat flux is determined from the 
bulk formulation. Interfacial salt flux can arise from bottom melting that balances 
heat flux and from a specified percolation velocity while conductive heat flux 
through the ice is set to zero (Eq. 3.15). From day 210 to day 215, the percolation 
velocity is set to LOxlO'*^ m s'', corresponding to a salt flux of 4.6x10'^ kg s''m'^ 
upward. This is the flux calculated from the change in salt content of the upper 
30 m. Before day 210, the percolation velocity is set to zero. The roughness 
length is set to 0.01 m. The bottom boundary conditions are the same as above. 

Figure 3.3 shows both modeled and observed salinity and temperature 
contours and elevation above freezing temperature at 5 m. The results are 
dominated by the enormous freshwater flux beginning on day 210. The fresh 
water is not distributed as deeply in the model as in the observations: on day 215 
the 30 psu contour is found at 11 m in the model and 17 m in the observations. 
The model output is adjusted to include the advective change in mixed layer 
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salinity from day 210 to day 215 so that the total salt content in both the model 
and the observed water column on day 215 are the same. The advective change is 
presumed uniform with depth and linear in time. The 30.5 psu contour is close to 
the observed, but the way in which fresh water is distributed above is quite 
different. The temperature contours show that relative to observations, the model 
ocean is colder at the surface and warmer at depth. Advective effects have not 
been used to correct the model temperature field, yet they probably are the reason 
for the variability in the extent of a warm core (-1.45 °C contour) during the 
period. The elevation of the mixed layer temperature above the freezing 
temperature drops steadily in the model. The model mixed layer became so fresh 
that the freezing temperature rose above the ambient temperature, implying the 
possibility of frazil ice formation. The supercooled mixed layer is visible from 
day 213 onwards in Fig. 3.3. The observed elevation above freezing temperature 
dropped from day 209 to day 212, then remained around 0.1 C° until day 215. 

There are four possible reasons that the observed mixed layer was warmer 
(relative to the freezing temperature) than was the simulated mixed layer. First of 
all, the model neglects all oceanic heating due to solar radiation. Of course, solar 
heating through leads and melt ponds was occurring, and even with no leads, 2- 
5% of incoming shortwave radiation enters the ocean through bare ice {Light et 
al, 2003]. Secondly, the SHEBA camp may have drifted into warmer water. It 
appears that the mixed layer and upper pycnocline were warming from day 212 to 
day 215 (Fig. 2.5), probably due to a combination of solar heating and advection. 
Third, the model did not mix the fresh water as deeply as the observations 
suggest, which means that warm water below the model pycnocline was not 
brought into the mixed layer and less thermal energy was available to maintain 
the elevation above freezing temperature. Underestimated vertical mixing plus 
the absence of the possibly advective warm core in the pycnocline hinder the 
model’s performance. Fourth, the model does not take into account frazil ice 
production such as that observed with the ROV in the upper 2 m during the initial 
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flushing of the lead. Such a release of latent heat and salt would to result greater 
temperatures and salinities deeper in the boundary layer. Unfortunately we do not 
have a quantitative observation of this phenomenon. 

Turbulence mast data from ocean city agree qualitatively but not in detail 
with one-dimensional model calculations. Figure 3.4 shows simulated friction 
speed and heat flux with ocean city TIC and AMTV observations. Salt fluxes 
from the TIC are not plotted since they are not reliable from 211-215. In Fig. 22, 
the TIC at 2.7 m very closely follows the one at 6.7 m (distance from interface) 
for both quantities, which is not generally true for the model. The model friction 
velocity at 6.7 m matches the observed for the most part. The model value of 
friction velocity at 2.7 m is nearly double the observed value on day 211, 
converging to it by day 213.5. The model compares more favorably to the AMTV 
values of friction velocity. The AMTV did not suffer from the same problems of 
biofouling as the TIC, nor is it vulnerable to possible biasing due to topographic 
features, as is the TIC. Biofouling and mechanical wear may have led to poor 
rotor performance, especially at 2.7 m. 

The observed heat fluxes are much different than the model calculation 
for days 212-214, where the observed heat fluxes were upward while the model 
heat fluxes are close to zero. The large upward TIC fluxes occur at about the 
same time that the surface stress momentarily increased and the AMTV observed 
even larger upward heat fluxes. It is possible that the inexplicably large upward 
heat fluxes at 5 m (AMTV) in localized areas halted the decline in the elevation 
above freezing temperature in the mixed layer on the whole, thus allowing the ice- 
ocean heat flux from day 212 to day 214 to remain strong. The other factors 
discussed above related to the low model heat content of Fig. 3.3 could also apply 
here. Ultimately, the large AMTV fluxes probably originated from the water 
heated in leads. Of course, the simulation does not include solar heating or the 
horizontal variability required for this mechanism. 
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The TIC data before day 219 certainly do not seem to be what is expected 
of a one-dimensional process driven by surface buoyancy flux. To investigate the 
effect of the massive influx of fresh water on the boundary layer structure and 
interaction with the ice, a 1-D, time-varying numerical model was used. The 
model was unable to correctly simulate the observations prior to day 219. A two- 
dimensional unsteady version of a three-dimensional model is now used to 
investigate the problem of the stable, stratified, non-homogeneous, unsteady 
boundary layer. First, performance in simulating days 219-220 is investigated, 
then the more extreme case of days 209-215. Perhaps it can explain the observed 
deeper, saltier mixed layer and abnormally strong upward fluxes on days 212-214. 

3.3 THREE-DIMENSIONAL MODEL 

A three-dimensional model has been developed with an eye towards 
quantifying the stability of the horizontally inhomogeneous boundary layer under 
melting sea ice. It is a non-hydrostatic, time-dependent, with the assumption of 
an incompressible, Boussinesq fluid. As discussed above, CTD observations at 
SHEBA from July and August show warm, fresh surface water mixing more 
deeply than predicted by a one-dimensional mixed layer model forced by 
observed winds and meltwater input. Autonomous vehicle observations around a 
lead during the same period show fluxes suggestive of a warm fresh lead plume 
underneath cooler, saltier seawater. It has been postulated that boundary layer 
shear in the presence of a horizontal density gradient would give rise to such 
unstable conditions [Crawford et al, 1999], increasing the vertical mixing beyond 
the one-dimensional prediction. It is the purpose of this section to explore this 
possibility. We also wish to explore the inherently two-dimensional, unsteady 
summer lead flushing process in its own right. 
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3.3.1 MODEL DESCRIPTION 

The model, a discretization of Eqs. 20-24, is based on the model of Smith 
et al. [2002] with modifications for summertime ice-ocean boundary layer study. 
Previous two-dimensional incarnations of the model are found in Smith and 
Morison [1993, 1998]. In Smith andMorison [1993], the upper ocean response to 
brine rejection in leads and movement of the ice was investigated. The two- 
dimensional hydrostatic model successfully simulated leads with a range of 
buoyancy flux and ice velocity, confirming the lead number scaling of Morison et 
al. [1992]. Smith and Morison [1998] showed that freely convecting wintertime 
leads show only minor departures from the hydrostatic simulations. It was also 
shown that two- and three-dimensional simulations allow for salty plumes to sink 
to the halocline relatively intact when ice velocity is low, actually shoaling the 
mixed layer when a one-dimensional model would deepen it. Higher-dimensional 
simulations can also allow for mechanical forcing at lead edges, which triggers 
convective events in their model. The most recent version of the model [Smith et 
al, 2002] is three-dimensional, non-hydrostatic with two options for 
parameterization of turbulent mixing. The generation of geostrophic eddies 
through dynamic instability around winter Arctic leads was examined. The model 
of Smith et al. [2002] forms the basis for the model described here. 

The conservation equations are those of Eqs. 3.1-3.5. The model is solved 
in a reference frame fixed to the ice surface, and a free-stream velocity term in the 
y-momentum equation (-fUi) forces horizontal flow. Stress is written in terms of 
local gradients. 


T = 


(3.18) 


As in the one-dimensional model the eddy viscosity is obtained as the product of 
friction speed and mixing length [McPhee, 1994]. Horizontally averaged surface 
values and vertical profiles of stress and buoyancy flux are used to specify the 
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eddy viscosity as a function of depth and time only. In the surface layer, friction 
speed is set to a constant (equal to the observed surface value) and mixing length 
increases linearly. Below the surface layer, friction velocity is determined locally 
from the shear (Eq. 3.10). The friction velocity is then horizontally averaged 
(being careful not to include the values in ice) and multiplied with mixing length 
to form eddy viscosity. Both horizontal and vertical momentum flux terms use 
the same viscosity, although they have differing background values. The 
diffusivity of temperature and salinity are equal to the eddy viscosity following an 
assumption of low gradient Richardson number. 

The pressure equation is derived as described in Smith and Morison 
[1998]. By taking the divergence of the momentum equation and applying the 
continuity equation an elliptic equation for pressure is obtained. 


dRHS^ dRHS ORHS^ 

- ^ +^ +^ 

dx dy dz 


V RHS 


(3.19) 


Here RHS is the right hand side of the vector momentum equation when 
only 5u / dt and the pressure gradient terms are on the left hand side. 

Centered differences in space and leapfrog in time conventions are used 
on a staggered grid. The rate of change of velocity is calculated without the 
pressure term, then new salinities (densities) and temperatures are calculated. 
Now pressure can be solved and used to revise the velocity calculation. Boundary 
conditions are applied, and the solution is averaged with the previous solution 
before advancing to the next time step. 

Doubly periodic flow is specified at the vertical boundaries, and no flow is 
allowed through the top and bottom boundaries. The bottom boundary is a free 
slip boundary for horizontal velocity, and no flux for either temperature or 
salinity. The top boundary is a flat surface lead surrounded by sinusoidal ice 
topography. There is a no slip condition for horizontal velocity. To impose 
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realistic ice topography, temperature and velocity are held constant in the ice. 
Multiyear sea ice at SHEBA was isothermal at the freezing point during the 
summer [Perovich et al, 2003]. The salinity flux is zero at the ice-ocean 
interface—^no ice melting or forming. In the lead a radiative heat flux is 
distributed with an extinction depth of 4 m. The lead salt flux is also specified, 
and it enters in the top grid cell only. 

The 40-m wavelength of the ice topography is chosen from an estimate of 
the peak in the spectrum of ice draft gathered from the AMTV. The 2-m 
maximum thickness corresponds to a mean ice draft of about 1 meter, and typical 
ridges with 2 or 3 m of draft, as indicated by the AMTV data. The ice topography 
is greatly simplified by using only one sinusoidal component. 

Initial conditions are either a steady state boundary layer solution for the 
given free stream velocity (assuming neutral stratification and a flat upper 
surface) or uniform velocity relative to the motionless upper surface. 
Temperature and salinity are initially uniform at 0 °C and 31 psu. 

3.3.2 RESULTS AND DISCUSSION—DAY 220 

The two most significant three-dimensional model experiments are 
discussed: strong fresh water flushing during boundary layer spin up (around day 
207) and freshening of a quasi-steady boundary layer (day 219.9). Model 
parameters are summarized in Table 1. Note that the y-dimension is much smaller 
than the other dimensions, essentially eliminating the y-derivatives. Therefore, 
from this point forward, we refer to this as the 2-D unsteady model. 

Figure 3.5 presents the quasi-steady result 4 hours into the simulation of 
day 219.9 (in which time the deep water traveled about 1500 m). Ice topography 
is outlined. Salinity (color shading) shows that the fresh water entering in the top 
grid cell is mixed downward, while at the same time it is carried downstream. As 
expected, the salinity is lowest at the surface at the downstream edge of the lead. 
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Table 1. Parameters for two three-dimensional numerieal simulations. Spin up 
refers to a period such as before day 210 of the SHEBA experiment, 1998, when 
the fresh layer in the lead was shrinking. Quasi-steady refers to a period such as 
around day 219, 1998, when there was no fresh water layer in the SHEBA lead. 


Model Parameters 

Spin Up 

Ouasi-steadv 

Lead Width (m) 

80 

480 

Domain Size (x, y, z m) 

(240, 4.5, 10) 

(1440, 9, 22.5) 

Ui (m s‘') 

0.07 

0.15 

Lead Salt Flux (10'^ kgs'*) 

4.8 

2.0 

Lead Heat Flux (W m'^) 

200 

448 

dx, dy (m) 

0.5 

1.0 

dz (m) 

0.1 

0.25 

dt (s) 

0.5 

0.5 

Background (m^ s'*) 

1x10'^ 

IxlO'** 

2 1 

Background Ax,y (m s' ) 

1x10'^ 

1x10'^ 


Temperature (white contours) follows a similar pattern, with the warmest water at 
the surface at the downstream edge of the lead. Isotherms are sloped differently 
than isohalines since the temperature source is distributed exponentially with 
depth, while only the water at the surface is freshened. Velocity {u and w vectors) 
shows a vertical shear on average, although the topography modifies the flow to 
satisfy the continuity equation. The simulation is not strictly steady because as 
long as the lead fluxes are present, the boundary layer will become fresher and 
warmer, eventually building a strong, shallow pycnocline. However, it is 
approximately steady in the sense that the wind and fresh water input did not 
change significantly over several hours (also assumed in the one-dimensional 
model). 

Horizontal transects of model temperature, salinity, friction speed, heat 
and salt flux (Fig. 3.6) compare very well to AMTV observations from day 219.9 
(Figs. 2.12-2.14). Temperature at 5 m in the lead shows a nearly identical 
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horizontal gradient to the AMTV data from day 220: a rise of 0.01 K per 100 m. 
Temperature deereases on a slightly steeper slope downstream of the lead edge in 
both the model and observations. Model salinity in the lead shows the freshening 
rate with downstream distance is weaker than observed, although absolute salinity 
from the AMTV is biased by a decreasing trend, as discussed previously. Friction 
speed at 5 m increases from lead to ice: 0.007 m s'' to 0.009 m s'' in the model, 
consistent with the AMTV observations. Heat flux at 5 m reaches about 130 W 
m'^ downward in the last 200 m of the lead, comparable to the AMTV lead 
averages for that run (98 W m'^ and 144 W m'^). Model heat flux at 5 m becomes 
upwards (nearly 100 W m') immediately downstream of the lead edge and decays 
to about half that value over the next 400 m. The covariance estimates of heat 
flux for all of segments 3 and 4 do not pass the significance test, although when 
using 100 m segments, two of the four pass (40 W m'^and 97 W m'^). The AMTV 
under-ice heat flux average is similar to a model transect between 5 m and 10 m. 
It should be noted that the AMTV did not travel beyond about 300 m downstream 
of the lead, making this downstream comparison with the model difficult. Model 
salt flux at 5 m approaches 1.2x10' kg m' s' , which is close to the AMTV lead 
values of 1.7x10'^ kg m'^ s'' and 2.4x10'^ kg m'^ s''. The model salt flux decays 
by more than 50% in 100 m downstream. AMTV observations also suggested 
weaker salt fluxes under the ice. 

The two-dimensional unsteady model transects share many features with 
the two-dimensional steady model results, with some important differences. In 
better agreement with the AMTV observations, the water properties return to 
background values much more quickly downstream of the lead in the 2-D 
unsteady model than in the 2-D steady model. The 2-D unsteady model heat flux 
at 5 m becomes upward immediately downstream, rather than after 2 km as in the 
2-D steady model, giving rise to much larger heat flux divergence between 5 m 
and 10 m. The constant temperature boundary condition under the ice results in 
more oceanic heat loss (by eddy diffusion) than the bulk heat flux condition used 
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in the 2-D steady model. The 2-D unsteady model lost thermal energy at an 
average rate of 80 W m' to the ice during the run shown, while the 2-D steady 
model was around 50 W m'^. The melt rate implied by the unsteady model heat 
flux is 2.6 cm day ’ (a salt flux about 20% of that in the lead). Observations of 
ablation rate at the SHEBA lead edge support the 2-D unsteady model heat flux 
results. On day 218 five ice thickness gauges from 6 m to 51 m from the lead 
edge (adjacent to our hut and downstream of the lead) measured an average melt 
rate of 3.8 cm day'*, while the SHEBA camp average was about 1.3 cm day ’ 
[Perovich et al, 1999]. 

AMTV heat flux observations show mixed comparisons with the two 
models, but seem to support the 2-D unsteady model. Under the ice the AMTV 
did not observe strong upward (2-D unsteady) or downward (2-D steady) flux, 
rather a weak downward heat flux. However, as noted, the unsteady model result 
between 5 m and 10 m is similar to the 5 m AMTV observations. The role of the 
molecular sublayer in the observed boundary layer may be important in affecting 
depth distribution of fluxes [McPhee, 1992] and should be examined with further 
study. Finally, salt flux at 5 m decreases downstream of the lead in the 2-D 
unsteady model, but increases in the 2-D steady model. AMTV observations tend 
to support the unsteady model, since it did not measure large salt fluxes for much 
of the under-ice transect. 

In general, the two-dimensional unsteady model agrees closely with 
observations from day 220, despite less realistic ice boundary conditions than the 
two-dimensional steady model. Topography and advective terms that are not 
possible in the steady model result in more rapid vertical mixing. Near the start of 
the run, for example, when the fresh water is first introduced, weak, unstable 
density gradients are found under the ice. This spin-up effect will be examined in 
the next section. The mechanical forcing associated with an explicit lead edge 
also causes enhanced mixing. Disagreements between the AMTV results and the 
2-D unsteady model might be resolved by allowing ice melt in the model. 



80 


3.3.3 RESULTS AND DISCUSSION—DAY 209 

Compared to the quasi-steady experiment, the spin-up problem is 
characterized by larger fresh water flux, smaller ice velocity, and no initial 
velocity shear. See Table 1 for a complete specification of the run discussed here. 
The salt flux in the model lead is set approximately to the flux estimated from the 
rate of decay of the fresh layer in the lead as observed by the ROV early in the 
storm (day 206.5 to day 210.5). 

Figure 3.7 shows the simulation after one hour. (Beyond one hour the 
model solution became numerically unstable.) The ice topography is outlined and 
is of the same dimension as the quasi-steady run. At this time, the lead was 100 
m wide. Salinity (shaded) shows the presence of a layer just over 1 m thick in the 
lead. As the layer is flushed past the ice ridges, an unstable density gradient 
develops. The fresh plume mixes downward as it is swept downstream, but 
dense, salty water remains in the ice depressions above. Temperature (white 
contours) shows the warmest water is at the lead surface near the downstream 
edge as expected. Under the ice the fresh water plume has a maximum 
temperature in its core a few meters below the ridges. Velocity (blue arrows) 
indicates strong shear in the top 5 m, particularly below the ridge keels. 

Buoyancy frequency (Eq. 2.10) indicates regions of low stability under the 
ice (Fig. 3.8). The model criterion for mixed layer depth (used in the 
determination of mixing length) is a buoyancy frequency of 1.5x10'^ s'^ and is 
contoured in black. The -1x10’’ s'^ contour is in white. The regions of unstable 
density gradient are visible in the ice cavities just downstream of the lead inside 
the white contour. Indeed, “overrunning” of dense water over light water is 
simulated. However, overturning is not observed in the simulated velocity field. 
A turbulent Rayleigh number can be used to investigate the likelihood of buoyant 
instability in the model given an unstable salinity gradient, after Smith and 
Morison [1998]. 
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Ra = 


as ,4 
9 a — h 
dz 


(3.20) 


Although the criterion for instability is intended for quantifying the relative 
influences of molecular diffusion and turbulent heat transfer between two 
horizontal plates maintained at steady temperatures, it can provide insight into the 
current problem. In the turbulent heat transfer problem, when Ra is larger than 
about 10 , buoyant overturning convection is expected to take place rather than 
“diffusive” transfer. In this numerical model, the analog to molecular diffusion is 
turbulent diffusion used to parameterize unresolved eddies, and convective 
instabilities correspond to resolved overturns. The eddy diffusivity and viscosity 
are used rather than the traditional molecular values, and they are assumed equal. 
Here h is the distance over which the salinity gradient is found. Substituting 
values from the simulation in Fig. 3.7 (es/dz = 6x10' psu m' , h=\ m As = = 

SxlO'"^ m^ s'^) results in a value of 200, which is nearly large enough for 
instability. If the gradient were larger and/or the depth over which it occurred 
were larger, a resolved instability could occur. For example, if h=2 m, Ra= 3000. 
If the model could be run for a longer time, it is possible the salinity gradient 
could increase to a critical value at which point we could expect the model to 
show overturning. The model domain would have to be extended sufficiently for 
the plume to descend to the critical value before the next lead was reached. From 
Fig. 3.7, the plume appears to descend at a rate of about 1 m per 80 m. In 
addition, if the eddy viscosities were lower by a factor of two, overturns would be 
expected in the model. 

Although the model simulation in Fig. 3.7 only lasted one hour, it is 
worthwhile to compare the snapshot in time with observations and 1 -D (no lead) 
simulations from days 207-215 at SHEBA. Thermal energy is found at greater 
depths in the 2-D unsteady model. The salt and heat content in each horizontal 
layer of the 2-D unsteady model are compared with an equivalent 1-D simulation 
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in Fig. 3.9. The 1-D model is run with the 2-D model’s horizontally averaged salt 
flux (1.6x10'^ kg m'^ s’'), heat flux (67 W m'^), and surface friction speed (0.005 
m s'^). An extinction depth of 4 m is used to distribute the heat flux. The run is 
repeated with a freezing point surface boundary condition replacing the zero-heat 
flux condition. The temperature and salinity profiles from these two runs are 
averaged so that the integrated heat content agrees with the 2-D model. 
Specifically, the profiles from the constant surface temperature run are weighted 
twice as heavily as the zero heat flux run since the 2-D model has twice as much 
ice boundary as lead boundary. The difference in heat content between the 2-D 
model and the net input due to radiation can be used to calculate the average heat 
flux to the ice in the 2-D model: 7 W m'^. The 2-D model does not allow this 
heat flux to melt ice and cause a salt flux. This would affect the solution only 
slightly since the implied melt rate is small compared to the specified salt flux in 
the lead. To make a fair comparison, melting was not allowed in the 1-D 
simulations here—^the salt inventory is in agreement as well. 

The vertical distributions of thermal energy and salt differ. Maximum 
heat content is found at 4 m in the 2-D model, but at 1.5 m in the 1-D average. 
Reduced salinity water is found as deep as 5 m in the 2-D model, compared to 4 
m in the 1-D model. The heat and salt distributions in the 2-D case reflect deeper 
mixing of the surface water. The depths influenced by the lead are still quite 
shallow compared to observations, although a longer run could possibly allow the 
warm, fresh water to mix more deeply. For example, the 1-D result from Fig. 3.3 
shows the fresh water introduced on day 210 (traced by the 31 psu contour) finds 
its way to 11 m by day 211, and in the 1-D model in this section it only reaches 4 
m. Perhaps the 2-D model would show a similar or deeper fresh water signal if it 
could be run over several days. The 2-D spin up simulation also suggests a longer 
run would result in a more vertically uniform water column (compared to 1-D 
simulation in Fig. 3.3), particularly if the regions of unstable density gradient 
began to overturn. Although neither the 1-D nor the 2-D unsteady model seems 



83 


to predict the relatively deep observed mixed layer, the 2-D unsteady model 
suggests that horizontal variability in heat and salt flux, combined with ice 
topography may contribute to enhanced vertical mixing. 



HF (W m"^) u* (m s'^) S (psu) T (deg C) 
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Figure 3.1. Two-dimensional steady simulation of SHEBA lead on day 220, 

1998. Horizontal transects at surface, 5 m, and 10 m. The lead is in the upper left 
and the ice is moving from right to left. 
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Distance upstream of lead edge (km) 


Figure 3.2. Density contour for same run as Fig. 3.1. Units are density minus 
1000 kg m'^. Black line indicates where horizontal density gradient changes sign. 
Velocity field is also shown (sub-sampled). 
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Figure 3.3. One-dimensional run for days 209-215, 1998. No leads or solar 
radiation are present. Fresh water surface flux of 4.6x10' kgm' s' begins on 
day 210. Salinity and temperature fields, and elevation above freezing at 5 m 
depth are shown. A uniform freshening rate of 0.4 psu over days 210-215 has 
been applied to the model salinity to account for observed advective changes. 
Observed contours are shown with dashed lines and bold numbers (from Fig. 2.5). 
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Figure 3.4. Friction speed and heat flux at two depths for the same run as Fig 3.3 
compared with TIC observations at 2.7 m and 6.7 m below the ice-ocean interface 
and AMTV observations. 
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Figure 3.5. Simulation after 4 hours of day 219.9 with a nonhydrostatic, three- 
dimensional model in 2-D unsteady mode. Salinity deviation is shaded, and 
temperature deviation is contoured. Velocity field is also shown; the velocity at 
the bottom is 0.15 m s '. Only every 10* velocity vector in the horizontal is 
shown. Ice topography (in white) is at a fixed temperature and position, and there 
is zero salt flux at the interface. Heat and salt fluxes at the lead surface are 448 
Wm'^ and 2x10'^ kg m'^ s’'. 
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Figure 3.6. Horizontal transects at surface, 5 m, and 10 m for same model run as 
Fig. 3.5. Lead is the 480 m segment beginning at 240 m. 
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Figure 3.7. Simulation of day 209 after one hour with nonhydrostatic, three- 
dimensional model in two-dimensional mode. Salinity deviation is shaded, 
temperature deviation is contoured. Velocity field is also shown, only one of 
every 10 vectors in the horizontal is shown. Flow is from left to right. 
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Figure 3.8. As in Fig. 3.7, but the square of buoyancy frequency is shaded. The 
black line is the 1.5x10' s' contour (the model mixed layer is defined by the 
average depth of this contour). The -1x10’’ s'^ contour is in white and separates 
positive and negative stability regions. 
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Figure 3.9. The salt eontent and thermal energy in each horizontal layer after 1 
hour from the 2-D model (blue) and the 1-D model (red). Each layer is 0.2 m 
thick. Two-dimensional results are from day 209 run in Fig. 3.7. The one¬ 
dimensional model is run with the horizontal average interfacial salt flux (1.6x10'^ 
kg m' s' ), heat flux (67 W m' ), and surface friction speed (0.005 m s' ). The run 
is repeated with a freezing-temperature surface boundary condition replacing the 
zero-heat flux condition. The solid line is the weighted average of the two (see 
text), equal in heat content to the blue line. 
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4. CONCLUSIONS 

This section discusses the possible ways in which this work might scale up 
to large-scale ice-ocean processes and numerical models. Our observations 
provide valuable insight on the seasonal evolution of the ice-ocean system. We 
find the snow and ice melt does not enter the ocean immediately as it forms nor 
does it uniformly cover the surface. Rather, meltwater is retained near the surface 
(which lowers the albedo). The water drains into leads where it can build up in a 
strongly stratified layer if ice speeds are low. The release of fresh water from 
leads can occur in a sudden, massive “flush” challenging notions of horizontal 
homogeneity in the boundary layer. A simple heat balance calculation below also 
suggests that horizontal homogeneity of ice-ocean heat flux is tenuous during 
periods of large lead fraction and floe perimeter. Large-scale models could use a 
lead number approach and keep track of lead fraction and floe perimeter to delay 
the input of fresh water appropriately. These parameters may also help in 
correcting the ice-ocean heat flux for horizontal variability scenarios. 

4.1 THE FATE OF THERMAL ENERGY IN THE UPPER OCEAN 

To put our observations in perspective it is useful to compute an energy 
budget for the upper ocean at SHEBA. Maykut and McPhee [1995] showed that 
solar heating of the ocean in the Beaufort Sea area accounted for a large portion 
of ice melt during AIDJEX in 1975. The calculation described in this section for 
SHEBA also shows that solar energy input to the upper ocean is nearly balanced 
by loss to the ice and ocean warming. However, it appears that the standard 
estimate of loss to the ice (Eq. 2.8) is not as effective in countering the solar heat 
input as suggested by the relatively small observed ocean warming through the 
summer. A cartoon (Fig. 4.1) shows the energy exchange in and out of the upper 
ocean reservoir that could result in a change in temperature: solar energy into the 
reservoir via three different pathways, sensible heat transfer to the ice (causing 
melting, conduction and/or change in ice temperature), lateral oceanic advection 
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of thermal energy, and sensible heat transfer at the reservoir base in the 
pycnocline. Here we assume the ice is isothermal at 0 °C and that the heat flux 
through the pycnocline is negligible. Our results are consistent with these 
assumptions. 

The fractions of area covered by ice, melt ponds, and leads changed 
significantly three times during June, July, and August, so the time series from 
Perovich [2002b] is approximated by the values in Table 2. Also shown is the 
mean downwelling shortwave and longwave radiation for each period from the 
SHEBA Project Office. 

Table 2. Area fraction of bare ice, melt ponds, and leads and mean downwelling 
shortwave and longwave radiation during the summer at SHEBA. 


Time 


Ice 

Pond 

Lead 

lo 

LWi 


(Date) 

fhi 


fid 

IW m'^) 


152-166 

1 June - 15 June 

0.97 

0.00 

0.03 

276 

277 

167-206 

16 June - 25 July 

0.75 

0.20 

0.05 

234 

295 

207-227 

26 July - 15 Aug 

0.60 

0.20 

0.20 

137 

299 

228-242 

16 Aug - 30 Aug 

0.78 

0.02 

0.20 

88 

294 


The albedo of leads, aid, is 0.066±0.007 on cloudy days \Pegau and 
Paulson, 2001]. The effect of wind waves on the lead albedo is pronounced at 
large zenith angles (during spring) in sunny conditions, but during June, July, and 
August the sky is typically overcast. 

The transmittances of bare ice and melt ponds as a function of wavelength 
were measured and reported by Light et al. [2003]. Before June 30 the bare ice 
was not strictly speaking bare. It had at least some measure of snow cover. No 
transmittance data is available for this transitional, snow-covered ice condition, 
and here we assume it is equal to the transmittance of bare ice. This is likely an 
over estimate, but the net effect is small. 
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Since the incident radiation was measured with a broadband detector, a 
spectrally-averaged transmittance for each type is formed: 




The integrals are over the wavelength range of strong solar radiation, roughly 350 
nm to 2500 nm. Since the spectral intensity was not measured, the Iqx function 
was taken from Grenfell and Perovich [1984], An average of clear and cloudy 
spectral intensities was used. Transmittance, Tx, was measured on several days 
and at a few sites during July and August, and an average over the sites was used 
for each surface type. Cases where a layer of snow (on ice) or skim ice (on melt 
ponds) was on the surface were not used. While this dramatically affects the 
transmittances, both were relatively rare from late June to mid-August. Minimum 
and maximum transmittances over all sites were used to calculate lower and upper 
bounds on the solar radiation input to the ocean and will be discussed later. 
Spectrally averaged transmittances were 0.15 and 0.03 for melt ponds and bare 
ice, respectively, with minimum values of 0.11 and 0.02 and maximum values of 
0.32 and 0.05. 

The net longwave flux is estimated from Eq. 4.2 [Ruffieux et al., 1995]. 

LW„^,=s(LW^-crT;) (4.2) 

Longwave radiation can only leave or enter the ocean through open water, as the 
ice is an excellent absorber of infrared radiation. The emissivity, s, is near one 
because of the very low lead albedo. We use a value of 0.95. The Stefan- 
Boltzmann constant, a, is 5.67x10'^ W m'^. The surface temperature of the 
lead is To, expressed in degrees Kelvin. We use the temperature in the lead layer 
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observed by Paulson and Pegau [2001] on day 203 (1.6 °C, a near maximum). 
The longwave emission from leads is then estimated to be 307 W m' . 

The area fraction of the surface types and the fraction of incident radiation 
that passes through each allow an estimate of the solar energy flux that reached 
the ocean. 

/, = +/,?"<■ +f,Al-au)}K (4.3) 

= (4.4) 


The amount of energy absorbed at the end of each time interval (Qnet) is shown in 
Fig. 4.2. By the end of August, the ocean absorbed 160 MJ m'^, roughly 60% 
from leads, 20% from melt ponds, and 20% from bare ice. Therefore, it is not 
appropriate to neglect the radiation transmitted through melt ponds and ice, at 
least during the SHEBA summer. The rate of solar heating of the ocean began to 
level off as the intensity of light diminished and as melt ponds froze over (mid- 
August). Neglecting the net longwave flux would result in a 5 MJ m' increase by 
day 227 and 9 MJ m'^ by day 242. 

The net radiative flux calculation is necessarily crude because of a number 
of factors. First, the transmittance is a difficult measurement to make, especially 
under melt ponds where edge effects are difficult to quantify [Light, personal 
communication]. Secondly, there are relatively few transmittance measurements. 
Certainly transmittance depends on a number of parameters that must be 
neglected for lack of data, among the most important are ice thickness and the 
type of surface. Thin snow layers or skim ice are very important, and even with 
more data, the detailed distribution in time and space of such surface coverage is 
very difficult to observe. Other sources of error include the slight mismatch in the 
wavelength range of the shortwave detector and transmittance data (285-2800 nm 
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vs. 350-2500 nm). However, sinee the intensity of radiation is weak below 350 
nm and above 2500 nm, the error in multiplying Iq by mean transmittance should 
be small. Also, other shortwave radiation measurements at SHEBA were 
collected and showed slight differences with the Project Office instrument. 
Biases of up to 5 W m' were observed [Persson et al, 2002], which would 
correspond to almost 40 MJ m' over June, July, and August. Persson et al. 
[2002] report that the Project Office instrument was biased low relative to other 
instruments, so we can regard our data as a lower bound on observed shortwave 
radiation. 

Net solar input is nearly balanced by basal ice melt (neglecting conduction 
and ice warming—both zero for isothermal ice held at the freezing temperature). 
The heat flux through the ice-ocean interface is estimated by Eq. 2.8, shown for 
part of the summer in Fig. 2.4. The net loss of energy through the interface 
during June, July, and August is calculated by integrating Eq. 2.8. 


Qice (0 = I pCpCuU^^ST dt' (4.5) 

h 


The temperature elevation above freezing is taken from the 5-m temperature and 
salinity observed with the ocean city CTD, which was far from the nearest lead 
(nearly 1 km). The friction velocity is calculated from wind and ice velocities 
using Eqs. 2.5-2.6. In Fig. 4.3, the energy lost to the ice is plotted along with 
solar input to the ocean (and its upper and lower bounds due to maximum and 
minimum transmittances). By the end of August, the ocean lost 130 MJ m' to the 
ice. The thermal energy lost to the ice estimated from average bottom ablation 
around SHEBA (green circles) is taken from Perovich et al. [2003] (their Fig. 
15b). They indicate less energy to melting early in the summer, and more energy 
to melting later in the summer compared to the bulk heat flux estimate. Overall, 
however, both results suggest ice melting was almost exclusively carried out by 
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locally heated surface waters, similar to the AIDJEX result of Maykut and 
McPhee [1995], 

In the simplest balanee, energy gained by the ocean but not lost to the iee 
would raise the ocean temperature. The ocean would have to inerease in thermal 
energy by 30 MJ m' by the end of August (cyan curve in Fig. 4.3), or about 10 
MJ m' using the lower bound of solar input (dashed cyan curve). The 
temperature of the oeean at four levels is plotted in Fig. 4.4 and does show a 
general warming near the surface for mueh of the summer. The thermal energy 
depends on the elevation above freezing (the upper set of eurves in Fig. 4.4), 
which is nearly identieal in shape to absolute temperature. The peaks in 
temperature early in the summer around 30-50 m seem to be unrelated to surface 
proeesses and are ignored here (perhaps a boundary current around Chukchi Cap). 
The temperature signal is very large and not evident at 5 m or 10 m. Conditions 
seem to return to normal on day 170. Then the temperatures above 30 m 
gradually elimb away from the freezing temperature until day 210. The strong 
winds after day 210 rapidly lowered temperatures before an interesting spike at all 
depths around day 217, when the camp drifted off Chukehi Cap. Later, the 
temperatures at 5 m and 10 m returned to near-freezing, while the 30 m 
temperature remained significantly above freezing. 

The ehange in thermal energy in the water eolumn above eaeh depth is 
calculated and plotted in Fig. 4.5. If one considers the ocean to only exehange 
energy at the upper surfaee, the curves in panel (a) are valid. The ehange in 
energy sinee June 1 (t = 0) in the water column above a depth, z, is given by 

0 

= \pC^[ST{t,z')-ST{0,z')-\dz' (4.6) 


If one assumes an advective heat flux into the upper ocean from adjacent regions 
can be represented by the change in heat eontent of the base of the mixed layer. 
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then this can be subtracted to obtain the change in internal energy due to surface 
processes alone. 

it, Z) = a, it, Z) - pC^\z\[ST{t-H) - dT{t)-H)-\ (4.7) 

The depth, H, must be below the level of influence of the surface, but small 
enough so that changes there are correlated to changes near the surface. We 
assume in Fig. 4.5b that the change in temperature at 30 m would have happened 
throughout the water column in the absence of surface processes. The advection 
correction does not make any difference between day 170 and day 210, as 
expected from the flat line temperature at 50 m in Fig. 4.4. After day 210, the 
energy in the upper layers fall in both formulations. Without correcting for 
advection, the thermal energy then increases above all depths (especially in the 
upper 30 and 50 m) from day 211 to day 220. Subtracting out the estimate of the 
advective change implies the water column is lower in thermal energy, i.e., the 
near-surface layers stayed the same temperature while the temperature rose at 50 
m. Presumably, the warmer water advected in was cooled by the ice, especially in 
the top 10 m. 

It is clear from Fig. 4.4 that some features observed at 50 m, and even 30 
m are not reflected at 5 or 10 m. The most obvious is from day 152-170, but there 
are other smaller events that make the advection correction unreliable and noisy 
(day 215, day 221, and 227). The assumption of homogeneous vertical change in 
the absence of surface processes does not always hold. The event around day 219 
seems to have been felt throughout the water column, as well as the general trend 
of warming from day 211 to day 220. It is unlikely that solar heating drove either 
signal since they were strongest at the deeper depths, while small but significant 
near the surface, which suggests these were advective phenomena. During this 
time, the SHEBA camp drifted from Chukchi Sea (500 m depth) over the shelf 
break to about 2000 m depth. It is likely that the complexity after day 210 in Fig. 



100 


4.5 is due to a combination of increased ice motion and vertical mixing as well as 
advection over a water column, which may contain a shelf slope boundary current 
whose properties are varying nonuniformly in the horizontal and even the vertical 
in the top 50 m. 

The cumulative change in thermal energy measured in the top 30 m of the 
ocean is plotted in Fig. 4.3 (red curve). The version “corrected” for advection is 
shown with the dashed line. In both cases, the energy in the water column from 
day 170 to day 210 increases at a rate slightly less than predicted from a balance 
of solar input and the estimated heat flux to the ice (cyan). An average oceanic 
heating rate of 2 W m' over the period would make up the difference. Using the 
minimum estimate of transmittance for melt ponds and ice would lead to 
agreement (cyan dotted) from day 170 to day 205. There are two possible reasons 
the ocean did not warm as much as predicted from the energy balance. First, as 
described, from day 170 to day 210, a fresh surface layer trapped solar energy in 
the lead. The surface layer was very warm (up to 2 °C), but this heating could not 
be measured by the ocean city CTD. Eight leads within 30 km of SHEBA had a 
warm, fresh layer on day 203 [Richter-Menge et al, 2001]. 

A second way to account for the lack of ocean warming is lateral melting 
of floe edges, which is not accounted for in the bulk heat flux formulation (Eq. 
2.8). Significant lateral melting was observed at two leads instrumented by 
Perovich et al. [2003]. In fact, Perovich et al. [2003] partitioned the total energy 
expended on melting into lateral melting (5%), basal melting (18%), and surface 
melting (77%) on July 20 (day 201). They also show an equivalent vertical heat 
flux for observed bottom and top melt rates to be 19 and 80 W m'^, respectively, 
implying an equivalent vertical heat flux for lateral melting of about 5 W m'^— 
about twice the value required to account for the difference between the red and 
cyan curves in Fig. 4.3. Using the balance based on minimum transmittance data 
implies agreement between the predicted and observed ocean warming, and 
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would mean no lateral melting or heating of upper leads took place from day 170 
to day 205, in disagreement with observations. 

The warm surface layers in leads were eroded beginning on day 207 (also 
see Richter-Menge et al, [2001]), but the mixed layer did not lose substantial 
thermal energy until the larger wind event on day 210. Then the observed ocean 
cooling is stronger than predicted, and the ocean remains cool rather than 
warming as the simple 3-way balance would suggest (cyan curve). A heat flux of 
15 W m'^ over the period would make up the difference. Either the solar input is 
overestimated (unlikely, see above), or the ice-ocean heat flux is underestimated. 
The latter is true when compared to the results of Perovich et al. [2003]. They 
show the energy used in bottom melting from June 1 to August 31 (day 152 to day 
242) was roughly 155 Ml m'^, compared to 130 Ml m'^ estimated from Eq. 2.8 
(Fig. 4.3). However on July 20 (day 201), the two estimates were much closer, 60 
and 50 MJ m' , indicating the problem with using Eq. 2.8 occurs during a time of 
high ice speed. However, given the wide applicability of Eq. 2.8, a more likely 
factor is a contemporaneous increase in floe perimeter and other boundary 
characteristics. 

The bulk heat flux formulation for bottom melt has been used successfully 
in the past [McPhee, 1992; Maykut and McPhee, 1995]. It assumes a horizontally 
representative mixed layer temperature and friction speed. However, Perovich 
and Elder [2002] observed large melt rates after day 213 near lead edges, and floe 
perimeter was observed to increase four-fold at the same time [Perovich et al, 
2002b]. While not previously considered to be important [Maykut and McPhee, 
1995], horizontal variation in mixed layer temperature over the scale of a single 
floe may be significant (see Fig. 2.12). The estimate of ice-ocean heat flux in this 
study could be low because after day 213 there was more floe edge, which 
experienced warmer temperatures due to adjacent leads. The AMTV observations 
(Fig. 2.12) show the temperature can vary by as much as 0.5 C° in 500 m. For 
this water to remain at the freezing point, it would have to freshen by about 1 psu, 
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which it does not. Therefore we can assume the freezing temperature behaves 
similarly to the in situ temperature. The horizontal variation in ffT is important, as 
is the horizontal variation in stress. 

The ice-ocean flux at ridge keels is likely to exceed the bulk estimate as 
well. A significant increase in roughness (say near ridge keels) may enhance the 
ice-ocean flux [Wettlaufer, 1991], and in fact the largest summer melt observed at 
SHEBA was an old ridge (average of 25 W m'^) [Perovich and Elder, 2002]. The 
relative importance of lead edges after day 210 seems to win out over keels 
however. While the effect of ridge keels is greater at the greater ice speeds after 
day 210, the observed increase of melt near lead edges was much larger then 
under keels \Perovich and Elder, 2002]. 

Generally, the results discussed here and in Perovich et al. [2003] are in 
good agreement. They both indicate the SHEBA melt season was characterized 
by extra solar insolation over previous years (at least 160 MJ m'^ versus 150 MJ 
m' at AIDJEX), important contributions from bare ice and melt pond 
transmittances (0.03 and 0.15 required for their heat budget to balance equal our 
estimate exactly), and the rapid increase in ice-ocean heat flux at the late-July 
storm (confirmed by ice melt and oceanic temperature drop). The melt rate was 
found to be high around lead edges, and the AMTV somewhat suggests the same. 

4.2 LARGE-SCALE IMPLICATIONS 

Eicken et al. [2002] point out that fresh meltwater is prevented from 
percolating through the ice into the ocean in the early summer by layers of 
superimposed ice, where fresh meltwater meets below-freezing ice or brine. Our 
observations of mixed layer salinity are consistent with this meltwater retention. 
Water that enters through flaws and cracks often forms false bottoms, which also 
prevents the release of fresh water into the ocean. Meltwater flowing into larger 
leads does not form false bottoms, and if the buoyancy flux is strong enough 
compared to the surface stress, the meltwater remains in a fresh surface layer. 
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The birth of this layer can be predicted from the Obukhov length, which can be 
calculated from large-scale parameters: ice-ocean stress and an estimate of 
buoyancy flux in leads. The latter flux is difficult to estimate because it is not 
simply the basal melt rate or the surface melt rate. It is strongly dependent on the 
permeability of the ice and the surface melt rate [Eicken et al, 2002]. As a first 
step in modeling such a process, one could assume that all of the surface melt runs 
into leads, i.e. impermeable ice. Then the surface melt would be multiplied by fbi 
fid^ to estimate the meltwater flux into leads. The surface melt could be held in a 
parameterized model surface reservoir until a critical Froude number was reached, 
at which point the lead flushing would begin. 

The effect of surface trapping of meltwater in a large-scale model would 
be to change the seasonal cycle of mixed layer depth. Currently, Arctic Ocean 
models such as the model of Zhang et al. [1998] use a bulk mixed layer with an 
entrainment velocity predicting the mixed layer deepening throughout the fall, 
winter, and spring. When the surface buoyancy flux is strongly stabilizing, the 
mixed layer depth jumps to the top grid point (say 10 m). Later in the summer or 
early fall, when the ice melt tapers off, the mixed layer depth begins to gradually 
increase again. Keeping the mixed layer deeper by withholding ice melt until 
later into the summer might modify this “saw tooth” pattern. The jump to 
minimum depth would be delayed until the flushing event. In this case, the model 
mixed layer would be fresher and shallower longer into the fall. Once the 
stratification from the flushing event was broken in early fall, heat stored before 
the flushing event would be available to melt ice or limit ice production. In other 
words, the sudden release of fresh water would put a tight lid on heat absorbed 
below the new, shallow mixed layer until a fall or winter storm could deepen the 
mixed layer. If the fresh water gradually came into the mixed layer, the ocean 
would not be capped off and the mixed layer could provide thermal energy from 
early in the summer onwards. 
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Currently it appears from AIDJEX (1975) and SHEBA (1997-1998) that 
the heat deeper in the mixed layer did not reach the surface until later in the fall or 
winter. We calculate roughly 5-10 MJ m'^ remain in the top 30 m from SHEBA 
observations, and Maykut and McPhee [1995] calculate 15-20 MJ m'^ for 
AIDJEX, mostly between 30 and 50 m. In effect, storage in the upper pycnocline 
weakens the ice-albedo feedback, since incoming solar radiation cannot reach the 
ice surface until later in the year when solar radiation is weaker. It would slow 
ice production somewhat. Large-scale models need to accurately simulate the 
fresh water flux in order to accurately simulate the ice-albedo feedback. Surface 
meltwater retention does two things: lowers the albedo early in the summer 
through melt pond formation, and it results in a strong fresh water flux that “seals 
off’ thermal energy in the pycnocline for future seasons. 

4.3 FUTURE WORK 

Small-scale, summertime, ice-ocean processes were observed in more 
detail than ever before in this study. A run-by-run analysis of the AMTV results 
may reveal even more information than was described here. The day-to-day 
conditions varied, and each run is a unique data set. The AMTV data set also 
contains very high-resolution ice draft data that has not been analyzed here. 

The small-scale processes have been simulated with one-dimensional and 
two-dimensional models. Both of the models used here were successful at 
reproducing many features of the horizontally inhomogeneous boundary layer. 
However, both models had difficulty with the lead flushing event on day 210 to 
day 215. In the 1-D model, an initial subsurface temperature maximum decayed 
away, and the mixed layer was colder, fresher and thinner than observed. The 2- 
D unsteady was more realistic, although it was only run for an hour due to 
numerical instability. It was unable to simulate the observed deep temperature 
maximum in the upper pycnocline, but it successfully simulated a shallower 
subsurface temperature maximum. It also suggested a new mixing mechanism 
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due to horizontal density gradient in the presence of vertical shear. This model 
could be improved by enhancing the numerical stability to allow longer run times. 
The role of the mixing mechanism could then be further investigated. 

A more sophisticated treatment of stratification and solar radiation may 
also improve the model. The model is initially neutrally stratified, so it is not 
possible to investigate the progression of the pre-existing mixed layer, and it is 
not possible to compare the modeled and observed heat content of the upper 
pycnocline. One way the ocean can warm in the pycnocline is by penetrative 
solar heating. The model uses a 4 m extinction depth for incoming shortwave 
radiation in the lead. However, the penetration of radiation in the visible is 
strongly dependent on chlorophyll [Kantha and Clayson, 2000]. In fact, incoming 
shortwave radiation can be attenuated to 10% of the surface value at a range of 
depths of 2-50 m, when chlorophyll ranges from 0.02 to 20 mg m' and the 
euphotic depth can vary from 10 to 100 m [Morel and Antoine, 1994]. At the 
SHEBA lead, chlorophyll was observed in the top 25-30 m, mostly due to two 
major phytoplankton blooms [Pegau, 2002]. Concentrations were as high as 1.5 
mg m' , which would correspond to a 90% reduction by 20 m and a 99% 
reduction by 40 m. This could allow on the order of 10% of the visible radiation 
to heat the pycnocline directly. A model pycnocline and a more sophisticated 
treatment of shortwave penetration may reveal under what circumstances a warm 
core is maintained, as throughout the storm of day 210. 

The quasi-steady period (day 220) was simulated very well, particularly 
by the 2-D unsteady model. It differed from AMTV observations in its prediction 
of heat flux under the ice downstream of the lead. However this may be due to a 
slight difference in the depth penetration of the boundary layer disturbance; the 
model results between 5 and 10 m agreed more closely with the AMTV 5 m 
results. The model agreed fairly well with ablation measurements at the lead 
edge. Overall the 2-D unsteady model results illustrate the importance of spatial 
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variability and the advective terms in exploring the boundary layer observations 
even under relatively steady foreing. 

Large-seale heat budget caleulations need to accurately represent the ice- 
ocean heat flux, which we assert is dependent on lead fraction and floe perimeter. 
While reasonably effective, the bulk exchange formulation does not agree exactly 
with the observed average melt rates in late summer. This is likely because actual 
rates vary greatly depending on ice type and location. Ridges and lead edges 
thinned the most, especially after the lead fraction and ice-ocean stress increased 
[Perovich et al, 2003]. The horizontal variability in ice-ocean heat flux needs to 
be incorporated somehow into any large-scale estimate. Horizontal variation in 
heat content and in ice-ocean stress are important factors to consider. 
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Figure 4.1. Illustration of the heat balance of the upper ocean under sea ice. The 
central box is the rate of change of thermal energy inside the dotted lines. This 
could be nonzero due to solar insolation through leads with a given albedo, or 
through melt ponds or ice with a given transmittance. The area fraction of each 
type is indicated by f Sensible heat flux at the ice-ocean interface and the 
pycnocline as well as advection of warm water also contribute. Net longwave 
radiation through leads is not shown; it is generally small unless lead fraction is 
high. 
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Figure 4.2. The accumulation of shortwave irradiance to the ocean through leads 
(red stars), melt ponds (green crosses), bare ice (blue circles), and the total 
(dashed line). Energy is contributed from melt ponds mainly during July, and 
from leads mostly from late July to early August (when lead fraction was highest). 
Including net longwave losses through leads results in the black line. See Table 2 
for area fractions and solar insolation values used to make this plot. 
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Figure 4.3. The major terms of the heat budget for the upper 30 m of the oeean at 
SHEBA. Net radiation input (blaek) is the same as Fig. 4.2, but now with dashed 
lines representing the maximum and minimum inputs based on the highest and 
lowest observed transmittances. Oeeanic loss (green line) is caleulated from the 
surface friction velocity and elevation above freezing at 5 m from the ocean city 
CTD. The difference in these two is the predicted ocean warming (cyan) and 
minimum warming (dashed cyan). The actual change in heat content is shown by 
the red (solid) with the dashed accounting for advective change in heat content. 
Heat used for bottom melting from balance measurements shown in green dots 
[Perovich et al, 2003]. 
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Figure 4.4. Temperature (bottom) and elevation above salinity-determined 
freezing temperature (top) at four different depths measured by the main ocean 
city CTD during June, July, and August. Data have been smoothed with a 1-day 
running mean. 
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Figure 4.5. Change in thermal energy of the ocean above four depths based on 
SHEBA ocean city CTD casts. In (a), the change is calculated by simply adding 
the integral from depth to the surface to the previous cast. In (b), the vertical 
integral of an imaginary 30 m column that changes uniformly at the rate observed 
at 30 m is subtracted from (a). Data have been smoothed with a 1-day running 
mean. 
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APPENDIX A 

Spectral Determination of Friction Speed 

Let the weighted spectral density of vertical velocity {(/>) be a function of 
turbulent kinetic energy production {P) and wavenumber {k). 


m^ks^{k)=f{p,k) 


(A.1) 


P = 


dU 

dz 


_ _ 3/ -1 


(A.2) 


In Eq. A.2, we assume shear times eddy viscosity equals the kinematic stress. 
The eddy viscosity equals friction speed times mixing length (A), so the shear 
turns out to be friction speed divided by mixing length. An empirical constant of 
proportionality approximately equal to 0.85 is c^. Using dimensional analysis to 
obtain the correct dimensions of </>, McPhee [2003] obtains the following 
relationship. 



r 


f—T' 

k 

V^max J 



(A.3) 


Therefore, the non-dimensional spectrum of vertical water velocity (weighted 
spectral density divided by the square of friction speed) has a universal shape and 
level in the inertial subrange. If the proportionality constant, Cy, is known, 
observed weighted spectra can be appropriately nondimensionalized to achieve 
the universal spectrum; i.e., friction speed is calculated to do just that. 

McPhee [2003] uses turbulent instrument clusters to calculate friction 
speed directly using Eq. 2.4 in 3-hour averages of 15-minute blocks of data from 
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November through March at ice station SHEBA. He then calculates the weighted 
spectrum of vertical water velocity for each 15-minute segment and smooths into 
logarithmic wavenumber bins, averaging over 3 hours, and finally 
nondimensionalizing by friction speed squared. After nondimensionalizing the 
wavenumber by kmax, all 387 3-hour nondimensional spectra collapse well to a 
single curve in the center of the expected scatter. Taking the logarithm of Eq. A.3 
and finding the value of <!> at a given y in the inertial subrange provides an 
estimate of Cy (0.48). McPhee [2003] states that Eq. A.3 neglects the turbulent 
kinetic energy sink of stabilizing buoyancy and assumes horizontally 
homogeneous conditions. Our study takes place in July and August, when 
stabilizing buoyancy and horizontal variability were often strong, but we will 
proceed anyway with this in mind. The effect of stabilizing buoyancy will be 
discussed briefly at the end of this appendix. 

Using turbulent instrument cluster data from the summer to form 
nondimensional spectra of vertical water velocity gives results similar to McPhee 
[2003]. However, the inertial subrange is found to be less well-defined; the 
spectra seem to fall off at a much lower logarithm of nondimensional 
wavenumber (0.2 vs. 0.7). In Fig. A.l are plotted nondimensional spectra of 
vertical water velocity from three different turbulence clusters: two levels at the 
interior of a floe (about 3.5 m and 7.5 m depths from days 210-215 and 219-220 
of 1998) and one from a lead edge (about 5 m in depth from days 219-220 of 
1998). Each spectrum represents a 3-hour (floe interior) or 1-hour average (lead 
edge). It appears that the deeper cluster and lead edge cluster give slightly 
different results than the upper cluster in the interior, perhaps due to a very 
shallow boundary layer and/or significant buoyancy flux at the cluster levels. It is 
also found that the clusters give slightly different constants of proportionality 
(0.33, 0.49, and 0.25 for interior shallow, interior deep, and lead edge, 
respectively), found by averaging the values from the spectra. At the shallow 
cluster in the interior we see excellent agreement for friction speed calculated 



121 


with the two methods: directly and only using the vertical water velocity spectra 
(Fig. A.2). Other clusters show less correspondence between the indirect and 
direct methods. The value of 0.48 was used for for all clusters. It seems the 
method is quite sensitive to the choice of Cy, since using the derived value for each 
respective cluster generally leads to a proportion between the two methods closer 
to one. The correlation coefficients between the two methods are unchanged 
however. 

Now Eq. A.3 can be used to calculate friction speed from the water 
velocity spectrum calculated for each 100 m segment of AMTV data. First, we 
calculate the weighted spectra as described (Fig. A.3). The inertial subrange, 
although small, seems to be centered on y* = 0.1. Note the double-humped 
structure in many of the spectra. All valid 100-m segments have been used, 
regardless of run depth. Runs resulting in a friction speed of less than 0.008 m s'* 
are in red. A three-dimensional view of the spectra illustrates the time 
dependence and shows higher spectral levels during the storm from days 210-215 
clearly (Fig. A.4). Once y* is chosen, we calculate the value of the 
nondimensional universal spectrum 0 ( 7 *) from Eq. A.3; we use cy=0.48, as in 
McPhee [2003], and we note that there is evidence to suggest smaller values. 
However, these values may be due to significant buoyancy flux or a violation of 
horizontal homogeneity. Using values of 0.25 or 0.33 introduces a uniform 
positive bias of 39% or 21%, respectively. We choose the value that matches the 
surface estimate of friction speed most closely, which also happens to be the value 
found by McPhee [2003]. Next we look up the value of the actual weighted 
spectrum log ^ at log k* = log kmwc +log y*. Finally, we use the definition of O to 
calculate the friction speed: 


(A.4) 


logM. = 0.5(log<zi, -logO(y.)) 
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A confidence interval is calculated from the confidence interval in the spectral 
estimate of (f) and amounts to a ±17% band. 

The procedure gives good agreement with surface friction velocity, despite 
the likely presence of stabilizing buoyancy flux and horizontal inhomogeneity, 
even for runs outside the boundary layer (Fig. 2.16). As long as the mixing length 
depends on the location of the maximum in the vertical water velocity spectrum, 
and there exists an inertial subrange, this method is valid in the boundary layer. 
As the buoyancy flux increases and mixing length decreases, the maximum in the 
vertical velocity spectrum moves to higher wavenumbers, and the inertial 
subrange becomes shorter. Eventually buoyancy flux could be so great as to play 
a part in the balance of the production of turbulence, and this method would fail. 

If we scale a new sink term in the turbulent kinetic energy equation, 
buoyancy flux, we see it has a similar form to the production term. 


Pb=-H^'p') = - 

P 


§ ^*0 
PkL^ 


(A.5) 


where we have simply applied the definition of the Obukhov length (Eq. 2.16). 
The Obukhov length plays a similar role in the buoyancy term to mixing length in 
the production term. As buoyancy flux becomes very large, both mixing length 
and Obukhov length approach zero, but their ratio approaches a constant. 
Therefore production and buoyancy flux destruction increase in tandem as 
Obukhov length decreases, which will continue until molecular scales are reached 
or horizontal homogeneity is broken by proximity to roughness elements. This 
method should work until that point because there still should be a region near the 
boundary, albeit small, where energy cascade depends on production and 
wavenumber (inertial sublayer or log layer). The key is making observations in 
that region that resolve the inertial subrange. 
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log 7 = log(k/k ) 

^ max' 


Figure A. 1. Vertical water velocity spectra from three turbulence instrument 
clusters deployed at SHEBA normalized by the covariance estimate of friction 
speed. The clusters are from depths of 3.5 m, 7.5 m, and 5.0 m. The former two 
were deployed at ocean city (mean ice draft of 1.4 m) from 210-215 and 219-220 
and the latter at the lead edge (“Sarah’s Lake”) on 205, 208, and 219-220. The 
slope of the inertial subrange is shown in green. By definition, all spectra have 
their maxima at log y=0. 












































124 



Figure A.2. Friction speed from spectral methods (ordinate) against friction speed 
directly from velocity component covariances (abscissa). The same three data 
sets from Fig. A.l are used. A proportionality constant (c^) of 0.48 has been 
assumed. Using the respective constant for each cluster leads to slopes closer to 
one, but the correlation coefficients remain the same. 
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Figure A.3. Weighted spectra of AMTV-derived vertical water velocity. Each 
100 m of data were used to calculate a spectrum, which has been fit with a 
polynomial to find the maximum. Only one of every five spectra is shown here. 
The friction velocity was calculated by the spectral levels at y=0.1, and the spectra 
have been color-coded according to the magnitude of the result (blue is above 
0.008 m s-'). 
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Dav of 1PP8 


Figure A.4. As in Fig. A.3, but all AMTV vertical water velocity spectra are 
shown in three dimensions. The color of the surface corresponds to its height 
(spectral density). Note the increase in spectral density in tandem with the 
interfacial friction speed (Fig. 2.16), which ultimately results in reasonable 
estimates of friction speed at the AMTV operating depths. 
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